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Abstract 


Electrophoretic and immunological analyses of genetic variation within 
and among 135 populations from throughout the range of the Plethodon 
glutinosus complex suggest that it is comprised of 16 groups that have 
achieved the species level of divergence. Problems associated with 
taxonomically recognizing forms that are genetically, but not always 
morphologically, differentiated are discussed. 


Part I 
Geographic Protein Variation 


Richard Highton 
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INTRODUCTION 


Modern methods of population genetic analysis make it possible to 
evaluate geographic variation in species and to analyze the patterns of 
their genetic substructuring. The slimy salamander, Plethodon glu- 
tinosus (Green, 1818), ranges from southwestern New England west to 
Missouri and south to Florida and Texas. A previous attempt to analyze 
patterns of geographic variation in the species (Highton, 1962b) was 
unsuccessful because of the paucity of variable morphological charac- 
ters. The present electrophoretic study of geographic variation in 22 
presumed genetic loci in 135 populations has provided the data needed 
to analyze variation in the species. The results indicate that the P. gluti- 
nosus complex, until recently considered a single taxonomic species, 
consists of 16 divergent units that are evolutionarily close to, or already 
have achieved, species status. 

Two pairs of sympatric, reproductively isolated species, so similar to 
P. glutinosus in standard taxonomic characters that they had been 
included in that taxon (Highton, 1962b), occur in two different regions. 
Widely distributed on the Cumberland Plateau is a form previously 
named P. kentucki (Mittleman, 1951). It is sympatric with P. glutinosus 
and its taxonomic validity was demonstrated by Highton and Mac- 
Gregor (1983) and Maha et al. (1983). A second sibling species, P. aure- 
olus Highton (1984), occurs in southeastern Tennessee and south- 
western North Carolina and is sympatric throughout its range with 
another member of the P. glutinosus complex, P. teyahalee Hairston 
(1950), which was also recognized as a full species by Highton (1984). 
He found that P. aureolus, P. teyahalee, and P. glutinosus are sympatric at a 
site in Polk County, Tennessee, with apparent reproductive isolation 
among all three forms. Data for salamanders from the latter site are not 
included in this study of 135 populations. The study includes data for 2 
sympatric, non-interbreeding species found at 6 of the 129 localities: 
P. glutinosus and P. kentucki at 4 sites, and P. aureolus and P. teyahalee at 2 
sites. 


4 Biochemical Evolution in the Slimy Salamanders 


This paper is one of a series of studies on electrophoretically detect- 
able genetic variation in salamanders of the genus Plethodon (Highton 
and Webster, 1976; Highton, 1977; Larson and Highton, 1978; Duncan 
and Highton, 1979; Highton and Larson, 1979; Highton and Mac- 
Gregor, 1983; Highton, 1984; and Wynn, 1986). 


MATERIALS AND METHODS 


An attempt was made to obtain 25 to 30 individuals from each of the 129 
sites (up to 5 hectares of continuous woodland per site) scattered 
throughout the range of the species (Fig. 1 and Appendix 1). It was not 
possible to obtain the full number from 8 of the localities, but at least 10 
salamanders were collected from all of them. Genetic variation is low in 
the populations with smaller sample sizes, so genetic distance estimates 
are probably not affected significantly. Two separate species were 
represented at 6 of the sites, therefore additional material was obtained 
to bring the sample sizes of both forms up to an effective number. Large 
sample sizes are needed to estimate accurately the geographic patterns 
of allele frequency variation in polymorphic loci. 

Populations were selected to include (1) sites in all of the physio- 
graphic provinces in which the complex is known to occur, (2) material 
from or near type localities of all species or subspecies presently 
recognized or synonymized under the name P. glutinosus, (3) several 
populations of each eastern geographic variant described by Highton 
(1970, 1972), and (4) sites where sufficient numbers of salamanders 
could be obtained for a satisfactory analysis of genetic variation. 
Electrophoretic data are available from over 200 additional populations. 
Although sometimes based on smaller sample sizes or fewer than 22 loci, 
or both, these data have been useful in better delineating geographic 
ranges of groups detected in populations from the 129 localities for 
which we have more complete data. 

Genic variation was analyzed in 22 presumed independent genetic loci 
using the protein-buffer combinations listed in Table 1. These are the 
same loci used by Highton and MacGregor (1983) and Highton (1984). 
The genetic. data are listed in Appendix 2, except for one general 
protein locus (Pt-3) that is monomorphic in all samples. 

Salamanders were brought to the laboratory alive and their blood was 
centrifuged to obtain plasma samples. Blood albumin was scored on 7 
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percent acrylamide disc gels and transferrin on 10 percent gels using 
the method of Davis (1964). Whole animals were homogenized and then 
centrifuged to obtain aqueous extracts of proteins, which were stored 
below —70°C until utilized for electrophoresis. Electrostarch and 
Sigma starch were employed using the methods of Brewer (1970), 
Selander et al. (1971), and Shaw and Prasad (1970). When two loci were 
associated with a protein, the more anodal locus was assigned the 
numeral “1” and the more cathodal was designated “2.” Allelomorphs, 
hereafter referred to as alleles, were designated alphabetically by order 
of electrophoretic mobility, with “a” being the most distant from the 
origin. 


Table 1. Buffers used in electrophoresis and proteins assayed on each buffer. 


Buffer Assay 


Tris Versene Borate | a—Glycerophosphate Dehydrogenase (a—Gpd), 
Glutamate Dehydrogenase (Gdh), General Protein 
(Pt-1), 6-Phosphoglutonate Dehydrogenase (6-Pgd), 
Phosphoglucose Isomerase (Pgi) 


Lithium Hydroxide _ Esterase (Est), Leucine Aminopeptidase (Lap), 
Peptidase (Pep), General Protein (Pt-2), 
Phosphoglucose Isomerase 


Poulik Esterase, Fumarase (Fum), Lactate Dehydrogenase 
(Ldh), General Protein (Pt-3) 


Tris-Citrate pH 6.7 Glutamic Oxaloacetic Transaminase (Got), Lactate 
Dehydrogenase, Malate Dehydrogenase (Mdh), 
Phosphoglucomutase (Pgm) 


Tris-Citrate pH 8.0 Glutamate Dehydrogenase, Glutamic Oxaloacetic 
Transaminase, Isocitrate Dehydrogenase (Idh), 
Leucine Aminopeptidase 


Tris-Glycine (disc) Albumin (Alb), Transferrin (Trf) 


Unless otherwise indicated, estimates of genetic similarity are ex- 
pressed as Nei’s J, the normalized identity of genes, and genetic distance 
is expressed by Nei’s standard genetic distance, D (Nei, 1972). These 
measures were calculated from a computer program prepared and 
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furnished by Nei and Roychoudhury (1974). Other measures of genetic 
similarity and distance were calculated by a modified version of the 
BIOSYS-1 program of Swofford and Selander (1981), which they 
provided. Estimates of genic heterozygosity (H), the mean proportion of 
loci heterozygous per individual, were estimated from allele frequen- 
cies. Phenograms and distance Wagner trees were obtained by using the 
BIOSYS-1 program and by the Taxan program made available by Rita R. 
Colwell of the University of Maryland. 

This study, comparing alleles in 22 loci (some highly variable) in 
almost 4000 animals, took 8 years to complete. Because several tech- 
nicians and students assisted in the laboratory work and several differ- 
ent starch lots were used during the course of the research, the chances 
for error may be higher than in previous studies. However, samples of 
the homogenates of all animals were stored at or below — 70°C for the 
entire study so that every allele in all populations was run side by side 
with standards on comparison gels. A few populations (indicated in 
Appendix 2) had rare alleles that could not be compared with all 
standards because all of the sample had been depleted in previous runs. 

Two samples that were originally to be included in the study were 
omitted when it was concluded that they may represent hybrid popula- 
tions. The data have been transferred to other studies of hybrid zones 
among the various groups of the P. jordani and P. glutinosus complexes. 

For many of the populations the sample sizes for the two blood protein 
loci (albumin and transferrin) are much larger than those for the 
remaining loci because data were available from an earlier unpublished 
study of geographic variation in these two loci. 

Exact localities and elevations for the collection sites are presented in 
Appendix | and are shown in Figure 1. The samples are numbered in 
the order of their clustering on the UPGMA phenogram of J values (Fig. 
2). This phenogram is the only clustering procedure that includes all 
135 populations. Only 115 samples were used for the other tree- 
building methods. 

As in most electrophoretic studies, some genetic variation is probably 
present that could not be resolved satisfactorily in all individuals. For 
example, it was not possible to separate the albumins of groups 3 and 6 
from the middle Atlantic states, although in a 1969 unpublished study, 
Virginia Maiorana had consistently been able to distinguish the albu- 
mins of these groups by using the method of Smithies (1959). Thus 
variation in allele d of the Alb locus is obviously more complex than 
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D POPULATION 


Fig. 2. UPGMA phenogram of Nei’s / values. Sample 92 is indicated as being a 
member of Group 13 by the symbol “1I3A.’ Populations are numbered as in 


Figure 1 and Appendix 1. 
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indicated here. Other similar examples are probably alleles ¢ of the Est 
locus, b of the Gdh locus, b of the Idh-1 locus, b of the Mdh-1 locus, d of 
the Pep locus, c of the Pt-2 locus, and j of the Trf locus. Allele jj of the 
Est locus and allele h of the Pt-2 locus are null alleles. Allele e of the Pt-2 
locus is not present in any of the 135 populations of this study but is 
common in related species (Table 4). 

Additional specimens from most localities are preserved in the writer’s 
collection and will be deposited in the National Museum of Natural 
History (USNM) collection. 


RESULTS AND DISCUSSION 


Grouping of Populations 


Populations previously recognized as a single species, P. glutinosus 
(Appendix 2 and Fig. 2), have been found to have substantial genetic 
differentiation. Three different species of the complex (P kentuck, 
P. aureolus, and P. teyahalee) are sympatric with P. glutinosus and there is 
little or no evidence of hybridization among them (Highton and 
MacGregor, 1983; Highton, 1984). 

A number of populations are very closely related genetically with D 
values under 0.025 (IJ > 0.975, Fig. 3). Many of these populations are 
from the northern parts of the range and are probably derived from 
different source populations that dispersed northward when the last 
glaciers receded, as discussed below. Several populations from the Gulf 
Coastal Plain physiographic province just east of the Mississippi River 
are also closely related to each other. Most other adjacent populations 
have Nei genetic distances higher than 0.025. 

Arranging the populations into genetically related groups facilitates a 
discussion of their relationships. A number of divergent groups of 
populations are genetically quite homogeneous internally and are 
geographically cohesive (Fig. 4). 

Some major differences in the arrangement of groups are seen in the 
UPGMA phenograms and distance Wagner trees which use different 
measures of genetic similarity and distance (see below). Apparently, 
there is no a priori method of determining which of the various tree- 
building methods would provide a dendrogram that best indicates the 
evolutionary relationships of populations from biochemical data (see 
Nei et al., 1983, for a discussion of the problem). However, the various 
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Fig. 3. Lines connect populations of the Plethodon glutinosus complex that have J = 0.975. 


Samples 129-131 (P. aureolus) and 132-135 (P. kentucki) are omitted. 
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methods may suggest phylogenetic reconstruction hypotheses. There is 
also considerable information on morphological variation in the P. glu- 
tinosus complex (Highton, 1962b, 1970, 1972) that helps to interpret the 
genetic data. 

Any level of D chosen to separate populations into taxonomic groups 
is arbitrary. The lower the value of D selected, the greater the number of 
groups, and vice versa. However, there seems to be a level of divergence 
(D = 0.15) that groups the samples into units that not only are geo- 
graphically continuous but also are not morphologically divergent 
within groups. If a higher mean D level is used for determining the 
group limits, populations that are geographically separated or mor- 
phologically divergent, or both, are sometimes included in the same 
group. If a lower value of D is used, several groups that appear 
morphologically similar are subdivided into units that are not geo- 
graphically contiguous. Using a D level of 0.15 after studying the results 
is, of course, subjective and others might favor different methods of 
grouping (or might not group populations by using genetic distance 
data). However, genetic data obtained by electrophoretic evaluation of 
20 to 30 protein loci has clarified our understanding of the relation- 
ships of populations in many plethodontid salamander species (see 
Wake, 1981, for a review of the literature). The discovery of the clock- 
like nature of substitution rates in the proteins evaluated in electrophore- 
tic studies (see Sarich, 1977; Wilson et al., 1977, 1987) provides a strong 
basis for using genetic data in the analysis of evolutionary relationships 
of populations. If the amount of time that populations have diverged is 
correlated with a build-up of genetic divergence, then genetic distances 
should provide better estimates of phylogenetic relationships than does 
morphology. Morphology is well known to have evolved at highly 
variable rates (Simpson, 1953; Wake, 1981; Larson, 1984). Electro- 
phoretic, immunological, and DNA hybridization data in Plethodon are 
in good agreement (Mizuno and MacGregor, 1974; Highton and Lar- 
son, 1979; Maxson and Maxson, 1979), and support the hypothesis that 
independent biochemical methods consistently reveal evolutionary 
relationships among taxa. Morphological studies have generally pro- 
vided less satisfactory results, and several investigators have suggested 
different arrangements of the species of Plethodon into species groups 
(Highton, 1962b). If divergence time is also correlated with the 
evolution of reproductive incompatibility (Mayr, 1942, 1963), then tax- 
onomic decisions involving allopatric populations at the species level 
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might more likely be correct if genetic distance data are available. 

Baverstock et al. (1977) and Thorpe (1982) also have argued in favor of 
using a threshold of about 15 percent genetic differentiation as the 
amount of genetic divergence necessary for the recognition of allopatric 
populations as species. Baverstock et al. (1977) suggest that if two 
allopatric populations possess fixed electrophoretic differences at 15 
percent or more of their loci, then it is probable that they belong to 
different biological species. (Two populations that have fixed differ- 
ences at 15 percent of their loci would have an J < 0.85 or a D> 0.16.) 
In asurvey of the literature (excluding birds), Thorpe (1982) found that 
97 percent of J values between species are below 0.85 and that within 
species 98 percent of J values exceed 0.85. 

In this study there are two borderline cases. On the UPGMA 
phenogram (Fig. 2), the mean D for the 130 comparisons between 
samples 93-102 and 103-115 is 0.146. Because these two groups 
represent morphologically different populations from isolated geo- 
graphic areas, they are regarded as separate groups (see below). ‘The 
mean D value of 20 comparisons between samples 11-20 and 21-22 is 
0.155. Because these populations are all morphologically similar and 
come from contiguous areas, they are regarded as belonging to the same 
group. 

Of 1,198 within-group comparisons, 55 (4.5%) have D values above 
0.15. Fourteen of these involve comparisons of samples 21 and 22 with 
other members of Group 3, and nine are comparisons of sample 92 with 
other members of Group 13. Of 7,847 between-group comparisons, 444 
(5.7%) have D values below 0.15. 

The allele frequency data in Appendix 2 show a variety of differences 
among the groups. In about half the cases there are fixed differences in 
at least one locus between parapatric and sympatric groups so that all 
individuals of one group are distinguishable from all those of the other. 
In the remaining cases, loci that differentiate groups have only allele 
frequency differences. For example, the Gdh locus is the highest 
contributor to the differentiation of the parapatric groups 5 and 6 
(mean D = 0.18). Allele a has a high frequency in all populations of 
Group 5 and is fixed in 3 of 5 samples, whereas allele b has a high 
frequency in Group 6 populations and is fixed in 38 of 40 samples. 
However, one sample from Georgia of Group 6 (72) has a 0.38 frequency 
of allele a and two Group 5 samples (31 and 32) have frequencies of allele 
b of 0.24 and 0.23, respectively. Thus over 99 percent of the sala- 
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manders of Group 6 are different at the Gdh locus from over 90 percent 
of those in Group 5, but three populations have a substantial frequency 
of alleles from the other group. In some loci, for example the Pt-2 locus, 
which usually distinguishes the parapatric groups 6 and 13 (mean D 
between the two groups, 0.24), there is a difference between all 41 
samples of Group 6 (fixed for allele a) and 10 of 12 samples of Group 13 
(polymorphic for alleles f and h), but 2 samples (126 and 127) of Group 
13 are fixed for allele a. A third example is more complex. The 
sympatric species P. teyahalee (Group 1) and P. aureolus (Group 15) have 
major frequency differences at six loci and the mean D between the two 
groups is 0.42, but always at least one population of one or both species 
has in low frequency an allele characteristic of the other species. This 
phenomenon might be due to ancient polymorphisms with the alternate 
rare alleles still occurring in only a few living descendants of the 
common ancestor, rare hybridization between the two species, or 
hybridization of P. aureolus with Unicoi Mountain P. jordani, which at 
some loci has a high frequency of alleles also found in P. teyahalee 
(Highton, 1984). 

An UPGMA phenogram (Sneath and Sokal, 1973) of Nei’s J values of 
all 135 samples (Fig. 2) indicates a high level of genetic divergence 
among groups, not only of the four species already recognized but also 
in the remaining populations still assigned to P. glutinosus. Using the 
level of D = 0.15 for grouping populations, 16 major groups are 
indicated in the phenogram. Because the UPGMA clustering algorithm 
assumes equal rates of molecular evolution on all lines, an assumption 
which may not necessarily be correct, the evolutionary relationships of 
these groups will be interpreted below, after discussion of the distance 
Wagner trees (Farris, 1972). 


Group 1 (populations 1-6). Samples represent P. teyahalee Hairston 
(1950), already recognized by Highton (1984) as a distinct species 
because it occurs sympatrically with Group 15 (P. aureolus) at 28 known 
localities and with Group 6 (P. glutinosus) at a single site in Polk County, 
Tennessee, with little or no indication of hybridization. Adult sala- 
manders of this form may be distinguished from those of other groups 
by their large size and very small white dorsal spots. There is overlap 
with variation in this pair of traits only with populations of Group 12, 
but the latter usually have much larger white dorsal spots. Salamanders 
in Group 1 sometimes have very small red spots on the legs. Mor 
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phologically, groups 1 and 12 are primarily distinguished by the size of 
the white dorsal spots, but genetically, they are quite different (see below 
on the possible origin of P. teyahalee through hybridization of Group 12 
and P jordani). P. teyahalee occurs in the Blue Ridge physiographic 
province of North Carolina and in immediately adjacent areas of 
Tennessee and South Carolina west of French Broad River. In South 
Carolina it also enters the Piedmont physiographic province. The 
distribution of this form is shown in Figure 8 and in Highton (1970: Fig. 
5; 1984). There is very little geographic genetic variation throughout 
the range of this group. 


Group 2 (populations 7-10). Samples of populations 7, 9, and 10 belong 
to a form from the mountains of northeastern Georgia that is unusual in 
often lacking the dorsal spots so characteristic of almost all other 
populations of the P. glutinosus complex (except those from Jasper 
County, South Carolina, which also lack the lateral spots). Salamanders 
in Group 2 usually have abundant lateral white or yellow spotting. The 
distribution of this phenotype was mapped in Highton (1970: Fig. 5). 
Surprisingly, a nearby sample of salamanders with abundant brassy 
dorsal spotting (population 8) is also a member of this group genetically. 
Although Hairston and Pope (1948) assigned some specimens of this 
form to their P. jordani rabunensis Pope and Hairston (1948), (see also 
Hairston, 1950), the salamanders of Group 2 from northeastern Geor- 
gia are genetically different from the topotypic Rabun Bald population 
that is correctly assigned to P. jordani (Peabody, 1978). Therefore no 
name is now available for this group. 


Group 3 (populations 11—22). Samples of populations 12—16, 18, and 19 
are from the Coastal Plain physiographic province of Virginia, North 
Carolina, and northeastern South Carolina. Population 17 is from the 
eastern Piedmont physiographic province of Virginia. The range of 
this group also extends into the Piedmont of South Carolina (popula- 
tion 11) and population 20 is from Habersham County in northeastern 
Georgia. The Virginia populations of this form were recognized as 
distinct by their small size and very small, brassy dorsal spots (Highton, 
1972) and were referred to as the “Coastal Plain type.” The North 
Carolina populations are similar in appearance. The northern popula- 
tions (13-19) are virtually identical to each other genetically and have 
little genetic variation. Most populations from South Carolina have 
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more abundant lateral white or yellow spotting than do those from 
North Carolina and Virginia. The southwestern populations have 
larger, brassy dorsal spots (Highton, 1962b). Samples of populations 21 
and 22 are differentiated genetically from the remaining members of 
the group (D = 0.155). The type locality of Plethodon glutinosus 
chlorobryonis Mittleman (1951), at or near locality 14, is within the range 
of this group. The name Salamandra cylindracea Harlan (1825) probably 
should be assigned to a form represented by the population near 
Camden, South Carolina (Highton, 1962b). If Harlan’s description was 
based on animals from south of Camden, he probably was describing 
this form, but populations north of Camden are of Group 12 (see 
below). 


Group 4 (populations 23-29). Samples of populations 23-27 from the 
Coastal Plain physiographic province of southern South Carolina and 
two others from adjacent coastal Georgia just west of the Savannah 
River (28 and 29) are members of this group. Populations 24 and 26 
represent the unspotted morphotype from Jasper County, South Caroli- 
na, described by Neill (1948). The same phenotype also occurs as a 
common to rare variant in other South Carolina populations. Spotted 
animals of Group 4 are similar in appearance to southwestern popula- 
tions of Group 3. The old name Salamandra variolata Gilliams (1818), 
with type locality “southern States,” is available for Group 4 if Schmidt's 
(1953) restriction of the type locality to the vicinity of Charleston, South 
Carolina, is accepted. However, Gilliams (1818) stated that the speci- 
mens on which the description was based were received from the 
“Florida Party.” If they actually came from Florida, the name would not 
apply to this group but rather to Group 11. According to Weiss and 
Ziegler’s (1931) account of Thomas Says trip to Florida, Say might have 
collected the type specimen in South Carolina, the Georgia Sea Islands, 
or along the St. Johns River in Florida, or even on his overland trip from 
Philadelphia to South Carolina. There is nothing in the original 
description to aid in assigning the name to any of the P. glutinosus 
complex groups. 


Group 5 (populations 30-34). Samples of this group are from east- 
central Georgia. Salamanders of population 33, from the Piedmont 
physiographic province, are large in size, but the other four samples are 
from the Coastal Plain physiographic province and are smaller in size. 
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They have brassy dorsal spots but often fewer than do most other 
spotted P. glutinosus. No name is available for this form. 


Group 6 (populations 35—74). This group occurs over a wide area from 
eastern Alabama and northwestern Georgia north to Illinois and east to 
West Virginia, western Virginia and Maryland, Pennsylvania, New 
Jersey, New York, and southwestern Connecticut. The distribution of 
this form in the middle Atlantic states (Pennsylvania, Maryland, Virgin- 
ia, and West Virginia) corresponds to that of the “brassy-spotted type” 
whose range was mapped by Highton (1972). Salamanders of most 
populations are large sized with abundant, large, brassy dorsal spots 
and have dark chins, compared to sympatric light-chinned, smaller 
P. kentucki (Highton and MacGregor, 1983). Populations from eastern 
Kentucky northeast to New York (47—70) are all very similar to one 
another genetically. Other western and southern populations show 
much more local genetic differentiation. The name Salamandra glu- 
tinosa Green (1818), with type locality Princeton, New Jersey, clearly 
refers to this form. The isolated record from southern New Hampshire 
(Highton, 1963) most probably belongs to this group. 


Group 7 (populations 75—76). Samples are from Round and Kiamichi 
mountains, Oklahoma, and represent a type whose genetic distinctness 
has already been noted by Duncan and Highton (1979). The sample 
from locality 8 of their study represents the same population as that of 
locality 76 of the present study (although different individuals are used 
here). Group 7 appears indistinguishable morphologically from other 
nearby populations except that in Group 7 the melanophore pigmenta- 
tion on the chin is often reduced compared to the pigmentation on 
adjacent populations of Group 13. Because of the lighter chins of 
Kiamichi Mountain P. glutinosus, Blair and Lindsay (1965) suggested 
that this group hybridizes with sympatric P. ouachitae, but Duncan and 
Highton (1979) could find no genetic evidence to support this hypothe- 
sis. No name is available for this form. 


Group 8 (populations 77-87). In the Coastal Plain physiographic 
province of western Kentucky and Tennessee, western Alabama, Missis- 
sippi, and the Florida parishes of Louisiana, there is a small-sized form 
represented by population samples 77 and 79-87. It has no distinguish- 
ing morphological features and often resembles Group 11 in appear- 


aa) 
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ance. Sample 78 is from the Appalachian Plateau physiographic 
province of northern Alabama and although adults attain a larger size 
than those of the Coastal Plain populations, it is genetically similar to 
them and is assigned to this group. No name is available for this form. 


Group 9 (populations 88—90). Samples from the Coastal Plain physio- 
graphic province of central Louisiana and southern Arkansas are small 
in size and have distinctive dorsal coloration: large dorsal spots with 
diffused brassy pigment. No name is available for this form. 


Group 10 (population 91). The sample from southeastern Oklahoma is 
closest genetically to Group 9. In size and color pattern, however, it is 
similar to other nearby populations of groups 7 and 13. No name is 
available for this form. 


[Although sample 92 clusters with 10 samples from the southeastern 
Coastal Plain physiographic province (Group 11) in the UPGMA pheno- 
gram, it is closer genetically to samples of populations 126 and 127 of 
Group 13 (D = 0.08—0.09) than to any of the samples of Group 11. 
Therefore, this population is regarded as a member of Group 13 rather 
than of Group 11 and is shown as Group 13A in Fig. 2.] 


Group 11 (populations 93-102). Populations 93-102 from the south- 
eastern Coastal Plain physiographic province of Florida, southern 
Georgia, and southern Alabama are characterized by their small size 
and large, brassy dorsal spots. The name Plethodon glutinosus grobmani 
Allen and Neill (1949) is available for this form if the older name 
Salamandra variolata Gilliams is properly applied to another form (see 
discussion under Group 4). 


Group 12 (populations 103-115). Samples of these populations are 
characterized by their large size and large white dorsal spots. This 
group ranges from the French Broad River in western North Carolina 
and eastern Tennessee northeast through most of the uplands of 
western Virginia. It also occurs in the panhandle of northeastern West 
Virginia and in parts of the Coastal Plain physiographic province of 
eastern Virginia. The distribution of this “white-spotted type” in 
Virginia and West Virginia was mapped in Highton (1972). The 
northern populations of this group (105-110) are genetically almost 
identical with one another. The name Salamandra cylindracea Harlan 
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(1825) is available for this form, if it is based on Group 12 animals from 
north of Camden, South Carolina (see discussion of this name under 


Group 3). 


Group 13 (populations 92, 116-127). This group includes all the 
remaining samples from west of the Mississippi River. As discussed 
above, sample 92 is also regarded as a member of this group. Group 13 
includes all the samples from Missouri and Texas and those from 
northeastern Oklahoma and northern and central Arkansas. No 
specimens from the two small areas where P. glutinosus occurs in the 
Coastal Plain physiographic province of eastern Texas are available, so 
the genetic relationships of these populations are unknown. Popula- 
tions of P. glutinosus from the Ozark Plateau and Ouachita Mountains 
are large, dark-chinned salamanders with large, brassy dorsal spots and 
closely resemble those of Group 6. The sample of population 127 from 
the Coastal Plain of southern Arkansas is probably similar, although too 
few animals have been collected there to accurately estimate the 
maximum size. In contrast, salamanders of the Texas populations 
appear to have a lighter ground color and smaller spots. Grobman 
(1944) recognized the latter as a subspecies, Plethodon glutinosus albagula, 
diagnosed by its lighter chin, but Highton (1962b) noted that some 
Texas populations have dark chins and that the eastern white-spotted 
populations (groups 1 and 12) also have light chins, as do groups 2, 15, 
and 16. Two of the three Texas samples (122 and 123) are of the dark- 
chinned form and the other (116) is of the light-chinned type. The 
dark- and light-chinned populations in Texas appear to be differenti- 
ated genetically (D = 0.09 and 0.13 between sample 116 and the other 
two). Indeed, the ten comparisons of Texas samples 122 and 123 with 
samples 117-121 (from Arkansas and Missouri) show that the former 
have lower genetic distances (D = 0.05) from the latter than they have 
from sample 116. Highton (1962b) suggested that the Edwards Plateau 
of Texas might have been invaded twice by P. glutinosus populations 
from the uplands of Oklahoma and Arkansas, and that the earlier 
invasion might have evolved into the light-chinned form in Texas. The 
later invasion by animals similar to the present populations of Arkansas 
and Oklahoma is now represented by the dark-chinned populations in 
Texas. The genetic data are consistent with this hypothesis. However, all 
these populations are similar enough genetically to be included in one 
group, for which the name P. albagula is available. 
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Group 14 (population 128). This population from the upper Coastal 
Plain physiographic province of Richmond County in eastern Georgia is 
distinct genetically from all other populations. It also has unusual 
pigmentation; the dorsal spots have much less brassy flecking than do 
spots of nearby populations of several other groups. No name is 
available for this form. 


Group 15 (populations 129-131). These samples represent a distinct 
species, P. aureolus Highton (1984). It is sympatric with P. teyahalee 
throughout its range and is also sympatric with P. glutinosus (Group 6) at 
a single known locality in Polk County, Tennessee, with little or no 
evidence of hybridization with either form. This small-sized species has 
abundant large, brassy dorsal spots. 


Group 16 (populations 132-135). This group represents a distinct 
species, P. kentucki Mittleman (1951), as demonstrated by Highton and 
MacGregor (1983). Its albumin has been shown to be immunologically 
distinct by Maha et al., (1983). This small species has a light chin and 
fewer, smaller, and less brassy dorsal spots than has sympatric P. glu- 
tinosus of Group 6. It occurs in the Cumberland Plateau of eastern 
Kentucky, northeastern Tennessee, southwestern Virginia, and western 
West Virginia. The four populations included here (from Harlan and 
Pike counties, Kentucky, and Wise County, Virginia) are different from 
some of the other samples evaluated genetically by Highton and 
MacGregor (1983). Three from their study (two from Summers County, 
West Virginia, and one from Dickenson County, Virginia) are especially 
divergent (D = 0.30) to the other samples of P. kentucki. Because I have 
not been able to obtain adequate samples of P. kentucki from these 
eastern localities, they are not included in this study. However, it 
appears from the few data available that P. kentucki may be divisible into 
two forms.at the level used here to separate groups. 

The distribution of groups 1 to 14 and the mean Nei genetic distances 
between parapatric groups are shown in Figure 5. 

Reconstruction of phylogenies from electrophoretic data is subject to 
many possible sources of error: laboratory errors in scoring gels, 
sampling errors in the collection of specimens and in the selection of loci 
to be evaluated, the partly stochastic nature of the evolutionary process, 
convergence of band mobilities of different genotypes, the possibility 
that natural selection might have a stabilizing or a balancing effect on 
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Fig. 5. Top. Numbers of the groups. Groups 15 (P. aureolus) and 16 (P. kentucki) 
are omitted. Bottom. Mean genetic distances (Nei’s D) between parapatric 


groups. 
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some loci, the nonmetricity of Nei coefficients, the likelihood of gene 
flow through hybridization between some groups, and the assumption 
of the UPGMA method of equal evolutionary rates along all lines. In 
light of these possibilities for error, it is remarkable that only one 
population (92) appears, on the basis of geography, misplaced on the 
UPGMA phenogram. 

Farris (1981) has argued strongly against the use of UPGMA pheno- 
grams of Nei’s J or D values for making phylogenetic reconstructions. 
Not only are the phenograms constructed on the assumption of equal 
rates of molecular evolution, but Nei distances are not metric and 
therefore may violate the triangle inequality. Farris suggests that there 
are conceptual problems with all of the available distance measures 
commonly used by those studying molecular evolution. On the other 
hand, Tateno et al. (1982) and Nei et al. (1983) used computer simula- 
tions to study the accuracies and efficiencies of different methods for 
constructing phylogenetic trees from molecular data under the assump- 
tion of equal rates of evolution. Although they found that no genetic 
distance measure or clustering method was consistently superior, in 
general the topologies obtained by the UPGMA method using Nei 
distances were recommended. Felsenstein (1984) has also defended the 
use of genetic distance measures in inferring phylogenies because they 
may represent the expected rather than the actual amount of change. 

Because of the very high cophenetic correlation coefficients obtained 
for phenetic clustering dendrograms in previous electrophoretic studies 
on Plethodon (Highton and Webster, 1976: r = 0.97; Larson and High- 
ton, 1978: r = 0.99; Duncan and Highton, 1979: r = 0.95; Highton 
and Larson, 1979: r = 0.97; and Highton and MacGregor, 1983: 
r = 0.97), it appears that the data provide consistent estimates of 
genetic relationships among groups. It is possible, however, that the 
high cophenetic correlation coefficients obtained in the above studies 
were partly due to the fact that some of the comparisons involved 
relatively divergent forms. A phylogenetic tree constructed by the 
method of Fitch and Margoliash (1967) for the data in Highton and 
Larson (1979) resulted in a tree very similar to that in the UPGMA 
phenogram. Summing all branch lengths from the base of the tree 
(Highton and Larson, 1979: Fig. 2) indicates very constant rates of 
molecular change along most of the lines leading to the 26 living species 
of Plethodon, a result not inconsistent with the molecular clock hypothe- 
sis of reasonably constant rates of molecular evolution. 
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On the other hand, the cophenetic correlation coefficients for the 
UPGMA phenograms (calculated for the 115 most different popula- 
tions, see below) of this study, 0.84 to 0.89 (Table 2), as well as those 
obtained in Highton’s (1984) study of several P. glutinosus complex 
populations, 0.89, are substantially lower than those of the above 
mentioned studies. Whether this low correlation is due to the fact that 
there are no highly divergent outgroups included in the studies of the 
P. glutinosus complex, or to other factors (such as hybridization between 
adjacent populations of closely related groups or unequal rates of 
molecular evolution), is unknown. 


Relationships Among Groups 


The BIOSYS-1 program (Swofford and Selander, 1981) has made it 
possible to calculate several different measures of genetic distance and 
similarity in addition to the commonly used measures of Nei (1972, 
1978). These measures include those of Rogers (1972) modified by 
Wright (1978), Edwards (1971; 1974), Prevosti (Wright, 1978), and two 
(arc distance and chord distance) of Cavalli-Sforza and Edwards (1967). 
Different clustering and tree-building algorithms are also available for 
use with various measures of genetic distance. Two were used: the 
phenetic UPGMA method (Sneath and Sokal, 1973) and the distance 
Wagner method of Farris (1972). The trees obtained by the latter 
method were rooted in two ways: (1) at the midpoint of the longest 
patristic distance, and (2) using the most divergent group (P. kentucki) as 
an outgroup. 

Unfortunately, it was not possible to use the BIOSYS-1 program for all 
135 samples because of a lack of computer capacity. By removing 20 
samples, however, the program could be used for the remaining 115 
populations. I chose to remove the 20 populations that are closest 
genetically to other samples in the study. Thus no population sample 
was removed that had a D of more than 0.005 to any other sample, so 
practically all of the geographic genetic variation present in the original 
135 populations still remains. Samples removed were 4, 13-15, 17, 18, 
56, 58, 60, 62, 63, 65, 90, 95, 106-110, and 113. The UPGMA 
phenogram of the J values of the 115 selected samples was compared to 
the complete phenogram with all 135 populations; the two are virtually 
identical except for a few minor changes in branching order for the most 
similar populations of Group 6. 
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The UPGMA phenograms of all the different genetic similarity and 
distance measures are similar with only minor differences in branching 
order of closely related populations. In all the phenograms, Group 16 is 
the most different and groups 1, 14, and 15 almost always separate from 
all other groups near the base of the phenogram. Groups 2, 3, and 4 
always cluster together as do groups 5 and 6, groups 9 and 10, and 
groups 11 and 12 (with one exception: on the Edward’s distance 
phenogram, Group 12 clusters with Group | rather than with Group 11; it 
probably belongs with Group | [see below]). Groups 7 and 8 do not 
always cluster together but are sometimes associated with several other 
different groups before they group together. 

Although all the UPGMA phenograms using the different measures 
of genetic distance are similar to one another, the distance Wagner trees 
often show major differences in the relationships of groups. For 
example, Group 16 (P. kentuckz) is the most divergent group on all 11 
phenograms, but on the five trees rooted at the midpoint of the largest 
patristic distance, Group 16 clusters most closely with Group 14 on four 
trees and with Group 15 on one tree. P. kentucki and its sister group 
cluster most closely to groups 2, 3, and 4 on two trees, Group 5 on one 
tree, Group 12 on one tree, and Group 13 on one tree. 

P. aureolus (Group 15) and P. teyahalee (Group 1) are either the second 
or third most different groups in 10 of 11 UPGMA phenograms. 
However, on the distance Wagner trees rooted at the midpoint of the 
largest patristic distance, P. aureolus clusters most closely to groups 1 
and 12 on one tree, Group 7 on the second, Group 16 on the third, 
groups I, 12, and 13 on the fourth, and groups 8, 9, and 10 on the fifth 
tree! On the other hand, the relationship of P. teyahalee is the same on all 
the trees; it is the sister group of population 115 of Group 12. This 
relationship is significant in understanding the genetic relationships of 
groups | and 12 as discussed below. Because of the many differences 
among the various trees and phenograms, it is not possible to suggest a 
phylogeny of the groups. Although P. kentucki and P. aureolus are the 
most differentiated forms and either one or both probably differenti- 
ated from the proto-glutinosus stock earlier than the other groups, so 
many different trees might be expected if most of the remaining groups 
differentiated from one another at approximately the same time. 
Under this hypothesis, groups | to 14 might have begun diverging from 
each other at a time when climatic conditions were conducive to the 
subdivision of the species into many geographic isolates. This pattern 
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of subdivision exists in the closely related species, P. jordani, whose 
range is now subdivided into 22 geographic isolates (Highton, 1970, 
1972). That some groups are slightly more closely related to others 
might be partly the result of hybridization and exchange of genes 
between parapatric groups and partly due to chance (in the case of 
allopatric forms). 

Although some parapatric groups are apparently more closely related 
than most other groups (e.g., groups 2-4, groups 5 and 6, and groups | 
and 12), probably most have been separated for a long time, perhaps 2— 
6 million years, based on the molecular clock estimates of Nei distances 
(Maxson and Maxson, 1979). Some of the groups may have had smaller 
ranges then they do now and might have been isolated from one another 
by the cold and dry climates of the Pliocene and Pleistocene (Davis, 
1983; Delcourt, 1979; Whitehead, 1981; Wolfe, 1978), but it is not clear 
how all of the forms might have been isolated. Some groups may never 
have received gene flow from adjacent groups when separated by such 
major barriers as the Mississippi River (between groups 6 and 13, 8 and 
9, and 8 and 13). The present distribution of the groups does not favor 
the isolation-by-distance model of differentiation (Wright, 1943, 1969) 
because there are usually abrupt changes in allele frequencies at the 
boundaries of groups rather than gradual clinal increases in genetic 
distance with geographic distance. Moreover, geographically distant 
groups are not consistently more different genetically from those near by. 

Thus I hypothesize that P. kentuck: diverged first, then P. aureolus, 
then slightly later all or most of the other 14 groups were isolated at 
approximately the same time. Subsequently, there probably has been 
occasional interbreeding between some of the parapatric groups and 
some influence of this hybridization is still present in the descendants of 
the affected populations. Because of the large amount of genetic 
divergence that has already taken place, the effects of gene flow have 
been minimal in the recent past and will probably be even less in the 
future. This trend would probably have been true even without the 
human disturbance of habitats, but the destruction of so many of the 
woodlands to which these salamanders are adapted will make gene flow 
between adjacent groups even less frequent in the future. The / and D 
values of the comparisons of all the groups are shown in Table 3. 

Table 2 shows four goodness-of-fit statistics for the different pheno- 
grams and trees: Farris’ (1972) “F,” Prager and Wilson's (1976) “F,” Fitch 
and Margoliash’s (1967) “percent standard deviation,” and Sokal and 
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Rohlf’s (1962) cophenetic correlation coefficient (rcs). The lowest values 
of the first three coefficients and the highest rcs should be seen when 
the trees best fit the distance data. The use of Fitch and Margoliash 
(1967) trees was not considered because Farris et al. (1982) have found 
that they are less efficient than distance Wagner trees, and most workers 
have not been able to find the best tree from the large number gener- 
ated. Fitch-Margoliash trees would be even more of a problem with the 
P. glutinosus data set because of the large number of populations. 

The BIOSYS-1 program (Swofford, 1981) also provides optimized 
distance Wagner trees to make the goodness-of-fit statistics comparable 
to those of the other clustering methods. All of these trees show much 
less variance in branch lengths than the unoptimized distance Wagner 
trees. Two of the unoptimized distance Wagner trees with the best 
goodness-of-fit statistics (see Table 2) are shown in Figure 6. The 
goodness-of-fit statistics are not very different from those of the 
UPGMA phenograms. 

All the populations of Group 1 (P. teyahalee) are closely associated on 
all of the distance Wagner trees, but instead of branching off early as an 
independent line as they do on the UPGMA phenograms, on all the trees 
Group 1 clusters with population 115 of Group 12. This arrangement 
probably indicates the evolutionary relationships of these two groups 
rather than the relationship indicated by the UPGMA phenograms. 
Groups | and 12 are the only large-sized, white-spotted forms and their 
ranges are parapatric, separated by the French Broad River. This river 
also appears to limit the ranges of P. cinereus and P. serratus (Highton 
and Webster, 1976). Population 115 is from the mountains near the 
headwaters of this river and it is the member of Group 12 that is closest 
genetically to Group 1. Group 1 has a lower average D value (0.20) to 
Group 12 than to any other group (Table 3). A comparison of P. teyahalee 
with other southern Appalachian large Plethodon by Peabody (1978) 
shows it to be more closely related to P. jordani populations of the 
Nantahala and Unicoi mountains (Table 4, group D) than it is to any of 
the southern Appalachian types of the P. glutinosus complex. P. teyahalee 
has common alleles at two loci that are unique (Pgic and Trf a), two that 
are fixed for alleles that are common in Nantahala Mountain P. jordani 
(Alb-c and Pt-2 6), and two that are common in Group 12 (Est-e and Got-1 h). 

Salamanders of Group 1 often possess small red spots on their legs, a 
character much more abundant in Nantahala Mountain P. jordani 
(Highton, 1962b). Both the genetic data and the morphological data 
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suggest, therefore, that P. teyahalee might have had a hybrid origin as a 
result of interbreeding between Group 12, white-spotted P. glutinosus, 
and red-legged Nantahala Mountain P. jordani, At the present time, P. 
teyahalee hybridizes extensively with Nantahala and Unicoi Mountain P. 
jordani (Highton, 1970) at all known contacts of their ranges, although it 
rarely hybridizes with the other geographic populations of P. jordani 
with which it is widely sympatric. The hypothesis of a hybrid origin of 
P. teyahalee thus is consistent with the genetic data, and its close genetic 
relationship to P. jordani may explain why it hybridizes frequently only 
with the Nantahala and Unicoi mountain populations of that species. 
Because of its hybrid origin, P. teyahalee might be expected not only to 
share far more genes with these populations of P. jordani, but the two 
might also have fewer reproductive isolating mechanisms. Judging 
from the small amount of geographic genetic variation in P. teyahalee, its 
origin by hybridization may have been a rather recent event, or it may 
have expanded its range only recently. 

The length of the branch leading to P. teyahalee on all distance Wagner 
trees indicates more molecular evolution on that line than on the 
branches leading to populations of Group 12. This result would be 
expected if there had been an infusion of foreign (P._ jordani) genes into 
the ancestor of P. teyahalee. By this hypothesis of the origin of P. teya- 
halee, the difference in the apparent rates of molecular evolution as 
observed in the distance Wagner trees would be more the result of 
hybridization than of differential substitution rates. By assuming equal 
rates of evolution on all lines and by averaging J values between groups, 
the UPGMA phenogram fails to indicate this evolutionary interaction. 
On the other hand, some of the distance Wagner trees indicate close 
relationship between groups based on similarities between single popu- 
lations of each group (perhaps due to gene flow) or similarities due to 
slower than average rates of molecular change rather than real genetic 
relationships (see below). 

The four populations of Group 2 cluster together on all distance 
Wagner trees just as they do on the phenograms. On the distance 
Wagner trees they are indicated as closer to Group 3 than to Group 4, 
whereas on the phenograms groups 3 and 4 are sister groups. The mean 
Dis 0.21 between groups 2 and 3, 0.25 between groups 2 and 4, and 0.20 
between groups 3 and 4. If there is some gene flow between the 
parapatric groups, then it might be expected that the two allopatric 
groups (2 and 4) might be slightly more different genetically. The three 
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groups are clearly more closely related to each other than they are to any 
other group. 

The southern populations of Group 3 have a considerable amount of 
geographic variation. As noted above, the seven northernmost popula- 
tions (13-19) of this group are almost identical genetically. These same 
populations have very low variability. If this region was recently 
reoccupied by salamanders from the coastal plain of southern North 
Carolina (the retreat of the Wisconsin glaciation began only about 
15,000 years ago), then they would be expected to be genetically similar, 
particularly because of their low variability. The most differentiated 
populations of Group 3 are those from the northwestern part of its 
range (20-22). 

The seven South Carolina and Georgia coastal plain populations of 
Group 4 (23-29) show a considerable amount of local genetic differen- 
tiation. The most different populations are the two from Georgia (28 
and 29) from just west of the Savannah River. One of these (29) is so 
different that it clusters with other groups on 3 of 10 of the distance 
Wagner trees. All individuals of the two samples from the west side of 
the Savannah River are spotted dorsally and laterally, but all the known 
populations from just east of the river in Jasper County, South Carolina, 
are unspotted. 

The five populations of Group 5 always cluster together. They are 
closest genetically to Group 6 (D = 0.18) but are more different from 
other parapatric populations of groups 4 (D = 0.36), 11 (D = 0.22), 
and 14 (D = 0.25). 

Group 6 is the largest group and consists of 40 samples (35-74). 
Genetically, the northern 32 populations (39-71) are closely related and 
24 populations (47-70) are very similar. The latter group probably 
spread northward only recently after the retreat of the Wisconsin 
glaciation made the area again habitable. Even now, after thousands of 
years of warming of the region, it has barely been able to enter glaciated 
territory in the extreme northern part of its range. This limited 
movement is in sharp contrast to the great distance that its congener, 
P. cinereus, must have dispersed northward (Grobman, 1944). P. 
glutinosus is more active during warm weather than sympatric P. cinereus 
(Highton, 1972), and the latter species ranges only as far south as North 
Carolina, compared to the extensive southern distribution of the P. glu- 
tinosus complex as far as the Gulf coast. Although none of the peripheral 
samples show any evidence of gene flow with Group 13 (groups 6 and 13 
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appear to be completely isolated by the Mississippi River), there is some 
genetic evidence of gene flow at the periphery of the contacts of Group 6 
with other parapatric groups. Sample 71 shows decided genetic 
influence of Group 8 and indeed clusters with that group on all but one 
of the distance Wagner trees. Samples 72 and 73 show some genetic 
influence of Group 5 and cluster most closely with the latter on 5 of the 
10 distance Wagner trees. Sample 97 (of Group 11) from the Coastal 
Plain physiographic province of southern Alabama shows considerable 
genetic influence of Group 6. Indeed, sample 97 is so close to Group 6 
that only a small change in its mean D value to the latter group would 
place it with that group in the UPGMA phenogram of Nei’s J values. 
Group 6 also hybridizes with Group 12 (Highton, 1972, and un- 
published data) and with Group | (unpublished data) where their 
ranges contact, although groups | and 6 have apparently evolved repro- 
ductive isolation at one site in Polk County, Tennessee (Highton, 1984). 

Group 7 includes only two samples (75 and 76) from the Ouachita 
Mountains of Oklahoma. It is different from other nearby populations 
as noted by Duncan and Highton (1979). Genetically, Group 7 is closer 
to groups 8 (D = 0.17) and 11(D = 0.18), east of the Mississippi River, 
than it is to parapatric western populations of groups 9 (D = 0.24), 10 
(D = 0.23), and 13 (D = 0.25) (Table 3). However, on 6 of the 10 Wagner 
trees it clusters with some of the above western groups. It is not clear 
why the populations from Kiamichi (sample 75) and Round (sample 76) 
mountains (all part of the same continuous ridge) from southeastern 
Oklahoma are so highly differentiated (Group 10 is also quite distinct) 
when Group 13 is so widely distributed from Missouri to Texas. There 
must have been considerable isolation of populations in the Ouachita 
Mountain region since the Pliocene, as also indicated by the even greater 
amount of differentiation of the salamanders allied to P. owachitae in 
that area (Duncan and Highton, 1979). 

Group 8 includes 11 samples (77-87) from the western Gulf Coastal 
Plain physiographic province of western Alabama, Mississippi, south- 
eastern Louisiana, and western Tennessee, and all are rather closely 
associated on both the UPGMA phenograms and all the distance 
Wagner trees. It also occurs in western Kentucky. Its closest relatives on 
most of the distance Wagner trees are groups 9 and 10, but on the 
UPGMA phenograms Group 7 is closest. 

Group 9 includes three samples (88-90) from the north-central 
Louisiana and southern Arkansas Coastal Plain physiographic prov- 
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ince. On all distance Wagner trees and phenograms this group is closest 
to Group 10 (population 91). Groups 9 and 10 in turn are closest to 
Group 8 on all the distance Wagner trees and the UPGMA phenograms. 

Group 11 consists of 11 populations (samples 93-102) from the south- 
eastern Coastal Plain physiographic province of Florida, southern 
Georgia, and southern Alabama. The UPGMA phenograms show it to 
be closest to Group 12 (D = 0.15), but it is also close to Group 8 (D 

= 0.17), Group 7 (D = 0.18), Group 13 (D = 0.19), and Group 6 (D 

= (0.20). In appearance the salamanders of this group are more similar 
to groups 6, 7, 8, and 13 than to Group 12. Indeed, the large-sized, 
white-spotted members of Group 12 are among the most different 
animals in appearance to Group 11. An examination of the distance 
Wagner trees, however, sheds more light on their relationships. On some 
of the trees (e.g., the Edwards distance Wagner tree in Fig. 6), the 
populations of Group 11 are all located near the base of the tree. This 
location suggests a slower rate of molecular evolution in this group than 
in the other groups. Thus, if through sampling errors or slower rates of 
change in these 22 loci, Group 11 has changed less than the others, it 
would be expected to cluster on the UPGMA phenograms with the 
group(s) that had also changed less from the common ancestor, in this 
case Group 12 and population 92 of Group 13. Its closest relatives are 
probably the adjacent members of the parapatric groups 5, 6, and 8. 

The relationships of the 13 populations (103-115) of Group 12 have 
been discussed under groups | and 11. As in the northern populations 
of groups 3 and 6, the six northern samples of this group (105-110) are 
also genetically similar to each other. The remaining, more southern, 
populations have a considerable amount of local genetic differentiation. 

The 13 populations of Group 13 (92, 116-127) are closer genetically to 
those of Group 11 (D = 0.19) than to any other group. This finding is 
surprising in light of the geographic proximity of groups 7, 9, and 10 (D 

= 0.25, 0.32, and 0.29, respectively) but is probably explained by the 
apparently slow rate of molecular evolution in Group II. 

Group 14, a single highly differentiated population (128), is closest 
genetically to the geographically nearby groups 4, 11, and 5 (D = 0.23, 
0.24, and 0.25, respectively). Two of the distance Wagner trees place 
this group closest to Group 5, four put it closest to Group 16, and in the 
remaining four it is found near the base of the line leading to all groups 
other than Group 16. The close relationship of groups 14 and 16, 
indicated by distance Wagner trees when they are genetically so differ- 
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ent(D = 0.43), indicates a problem with some of the distance measures, 
the distance Wagner method, or both, for use in phylogenetic recon- 
struction of electrophoretic data. 

The three populations of Group 15 (P. aureolus) are very different 
from all others. Although they are closest genetically to Group 8 (D 

= (0.24), they are very different from these populations. In the various 
distance Wagner trees Group 15 clusters with groups 1, 7, 8, 9, 10, 12, 13, 
and 16! Considering the high amount of genetic distance that has been 
built up in P. aureolus from all other groups, it is surprising that in only | 
of 10 trees it appears as an outgroup to all the other groups except P. ken- 
tucki (the Edwards distance Wagner tree in Fig. 6). The latter arrangement 
is probably a better indication of its phylogenetic relationships. 

Group 16 (P. kentucki) is by far the most different group on the 
UPGMA phenograms. As in the case of P. aureolus, it clusters with 
several different groups on the various distance Wagner trees rooted at 
the midpoint of the greatest patristic distance. 


Pleistocene Effects on Interpopulation Relationships 


Highton and Webster (1976) found that all samples of Plethodon cinereus 
from the northern three-quarters of its range (those occupying glaciated 
territory) were almost identical genetically. This genetic similarity may 
indicate that all were derived from the most northern, closely related 
populations that survived the Wisconsin glaciation and were in the most 
favored position to disperse northward upon its retreat. Larson (1984) 
and Larson et al. (1984) have discussed gene flow in these populations of 
P. cinereus. Since the Wisconsin glaciers began their withdrawal only 
about 10,000-15,000 years ago, probably little molecular evolution 
would have occurred in northern populations and they all might still be 
very similar to their common ancestor. In the P. glutinosus complex, 
genetic uniformity is also apparent in the most northern populations 
within several of the groups. Davis (1983), Delcourt (1979), and 
Whitehead (1981) review changes in the eastern forest types during the 
last 30,000 years. 

After the last Pleistocene glacial maximum, Group 3, whose north- 
ernmost range probably was limited to the Coastal Plain physiographic 
province of South Carolina or southern North Carolina, probably 
moved northward in the coastal plain to southern Virginia. It appar- 
ently was able to cross the James and York rivers, but this group is not 
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known north of the Rappahannock River or on the Del-Mar-Va Penin- 
sula, across Chesapeake Bay, where the habitat appears ideal for this 
form that is adapted to the coastal plain. 

Group 12 probably moved northward from northern North Carolina 
or southern Virginia in the Piedmont, Blue Ridge, and the eastern 
Valley and Ridge physiographic provinces until it reached the Potomac 
River, which it apparently was unable to cross except in Washington 
County, Maryland (Highton, 1972). It has also extended its range into 
the coastal plain of Virginia in the “northern neck” (between the 
Potomac and Rappahannock rivers), as well as the area just south of the 
Rappahannock River. 

Group 6 probably moved northward from the Appalachian Plateau of 
southern Kentucky and spread to Ohio, Pennsylvania, Maryland, New 
Jersey, and New York. Because Group 12 was unable to cross the 
Potomac River (except in Washington County, Maryland), Group 6 was 
able to migrate eastward into the Blue Ridge and Piedmont provinces of 
Maryland and Pennsylvania, provinces that are occupied solely by 
Group 12 south of the Potomac River. Apparently Group 6, unlike 
Group 12, is not adapted to life in the coastal plain, for it appears to be 
absent from that province in Long Island, New Jersey, Delaware, and 
Maryland (Highton, 1962b, 1972). 

Populations of Group 6 in Illinois (41 and 42) and western Indiana (40 
and 44) are genetically closest to those in western Kentucky (39 and 43) 
and southeastern Tennessee (45). It appears that these western 
populations of Group 6 crossed the Ohio River (in one direction or the 
other) during the middle or late Pleistocene, subsequent to the genetic 
differentiation of some of the major subgroups of Group 6. 

Northern populations of Group 13 are not as similar genetically to 
each other as are those of the other three northern groups (3, 6, and 12), 
but the three Missouri samples (118—120) and the northernmost Arkan- 
sas sample (121) are among the most closely related members of that 
group. The latter four populations are also very similar to the dark- 
throated Texas populations (122 and 123). Whether this close relation- 
ship is the result of recent colonization of Texas by populations from the 
northern part of the range of Group 13, or colonization in the opposite 
direction at the end of the last glaciation by populations from Texas, 
cannot now be determined. 

The only other groups containing closely related populations are 
groups 1, 8, and 9. Group 1, as suggested above, may be of hybrid origin 
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and this event may have taken place rather recently, thus leaving little 
time for subsequent local differentiation. The populations of Group 8 
that are all very closely related (82—86, from Mississippi and south- 
eastern Louisiana) are unusual in that most other southeastern coastal 
plain populations of groups 4, 5, and 11] have a considerable amount of 
within-group local differentiation. If the rise in sea levels during 
interglacial periods of the Pleistocene were responsible for exterminat- 
ing populations in low-lying areas, then those populations that result 
from recolonizations might show less local differentiation than those 
from adjacent higher areas. However, it is not clear why this pattern is 
seen in Mississippi but not in Alabama, Florida, Georgia, and South 
Carolina. Perhaps the southern populations of Group 8 only recently 
colonized that area because it was previously uninhabitable for some 
unknown reason. 


Analysis of Substitutions within Loci 


Variation in rates of substitutions at the structural gene loci analyzed in 
this study is shown in Table 4. Of the 22 loci, 4 (a-Gpd, Ldh(h), Mdh-1, 
and Pt-3) have the same common (usually fixed) allele in all groups. 
Because the same alleles are also common or fixed in the other species of 
the P. glutinosus group, these loci are presumably changing at very slow 
rates in these salamanders. Five loci (Alb, Est, Got-1, 6-Pgd, and Trf) 
appear to be evolving at very rapid rates, in agreement with the 
conclusion of Sarich (1977) that plasma proteins and esterases are 
among the most rapidly evolving loci in vertebrates. The remaining 13 
loci appear to be evolving at intermediate rates and are useful in 
attempting to determine the ancestral allele when compared to the six 
species most closely related to P. glutinosus and the other two species of 
eastern large Plethodon (P. wehrlei and P. punctatus) of the P. wehrle: group 
used as an outgroup. One or more substitutions have occurred in some 
of the intermediate or fast evolving loci in all 16 groups, as seen in Table 
4. The data on the five related species were obtained from comparisons 
of P. glutinosus alleles with the samples of P. jordani used by Peabody 
(1978), and for other species by Duncan and Highton (1979), and Wynn, 
Highton and Jacobs (1988, and unpublished data). Highton and 
MacGregor (1983) and Duncan and Highton (1979) found that P. yonah- 
lossee is genetically the most different member of the species group. 
The other five species (P. jordant, P. fourchensis, P. ouachitae, P. caddoensis, 
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and P. petraeus) also differ from all of the groups of the P. glutinosus 
complex. 

Which alleles of the five fast-evolving loci were present in the common 
ancestor of the group cannot be determined because of the large 
amount of variation now present. On the other hand, for the remaining 
13 variable loci apparently evolving at slower rates, the likely ancestral 
allele can be selected for the 16 groups of the P. glutinosus complex with 
the aid of the data on the related species used as outgroups (Table 4). 
For example, at the Gdh locus, allele b is very common or fixed in all 
forms with the exception of P. glutinosus complex Group 5. Thus it is 
very likely that the substitution of allele a occurred during the differen- 
tiation of Group 5, although two of the five samples of this group (31 and 
32, see Appendix 2) still have a low frequency of allele b. It is similarly 
possible to suggest the probable ancestral allele at some of the other loci 
evolving at intermediate rates. In the case of Pt-2, the probable 
ancestral allele (e) is not present in any of the 135 samples of the P. glu- 
tinosus complex. However, because allele e is not only present in the two 
outgroup species, but in seven forms of the P. glutinosus species group, it 
is likely that it is the ancestral type. 

For each P. glutinosus group species, the number of loci evolving at 
intermediate rates that have probably had electrophoretically detectable 
substitutions replace (or become commoner than) the ancestral allele is 
shown on the bottom line of Table 4. Since the time that all species had a 
common ancestor, the number of loci evolving at intermediate rates 
varies from one to five with groups 3, 5, 7, 10, and 15 having the most. 
Groups 11 and 12 have had fewer substitutions and are the same at all 
intermediate-evolving loci, and it is therefore not surprising that they 
have such a low D value because the contribution to D between groups 11 
and 12 is almost entirely due to the few fast-evolving loci. From the 
present data it is not possible to determine whether their similarity is 
because of close evolutionary relationship, a significantly slower rate of 
biochemical differentiation, or simply because the loci selected for 
analysis by chance happened to have had fewer substitutions. The first 
of these possibilities appears unlikely because of the large amount of 
morphological divergence between the two groups and their geographic 
separation, but the correct choice between the second and third 
possibilities probably could be made by further study of additional 
independent genetic loci. I think that the last explanation is probably 
correct. Attempts to do a cladistic analysis of substitutions, such 
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as attempted by Wake et al. (1978) and Patton and Avise (1983), 
have proven impossible because of considerable polymorphism at some 
loci. 

An attempt to find the relationships of the 26 forms of the P. glutinosus 
species group was made by comparing one sample of each form (usually 
the sample selected was from at or near the type locality) with the two 
species of the P. wehrlei group. UPGMA phenograms and distance 
Wagner trees all indicate that the 16 species of the P. glutinosus complex 
are not a monophyletic group, but instead cluster variously with other 
species of the species group. Only one pair of species of the group, P. cad- 
doensis and P. ouachitae, cluster together consistently, as do P. wehrlei and 
P. punctatus of the P. wehrlei group. This result is interpreted to indicate 
that the other species also diverged from the same common ancestor at 
approximately the same time as the 16 members of the P. glutinosus 
complex. 


Variability of Populations 


The average heterozygosity of each population sample is shown in 
Appendix 2 and indicated geographically in Figure 7. The means and 
ranges for each group are shown in Table 5. 


Table 5. Heterozygosity of groups. 


Group Species Population Mean Range 
1 P. teyahalee 1-6 0.05 0.02—0.09 
2 P. chattahoochee 7-10 0.19 0.14—0.23 
3 P. chlorobryonis 11-22 0.04 0-0.13 
4 P. variolatus 23-29 0.07 0.04—0.10 
5 P. ocmulgee 30-34 0.06 0.04—0.08 
6 P. glutinosus 35-74 0.07 0.01—0.18 
7 P. kiamichi 75-76 0.11 0.08—0.15 
8 P. mississippi 77-87 0.06 0-0.13 
9 P. kisatchie 88-90 0.002 0—0.004 
10 P. sequoyah 9] 0.04 — 
11 P. grobmani 93-102 0.06 0.01—0.12 
12 P. cylindraceus 103-115 0.05 0.02—0.10 
13 P. albagula 92,116—127 0.06 0.01—0.14 
14 P. savannah 128 0.0 _ 
15 P. aureolus 129-131 0.12 0.08—0.15 
16 P. kentucki 132-135 0.13 0.11-0.15 
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All but two of the smallest groups appear to have a considerable 
amount of genetic variation. Group 2 has by far the highest amount of 
variation with a mean H = 0.19 and includes the two most variable 
samples with H = 0.23. Groups 7, 15, and 16 also have a large amount 
of variability with mean H ranging from 0.11 to 0.13. The four samples 
of groups 9 and 14 have extremely low genetic variability. The mean H 
values of the remaining 10 groups indicate a moderate amount of 
genetic variability (0.04—0.07). The unweighted mean H of all 135 
populations is 0.06. 

With a few exceptions, there does not appear to be a consistent pattern 
of high or low heterozygosity at the periphery or in the center of the 
ranges of the various groups. This pattern agrees with that in Drosophila 
summarized in Brussard (1984). The northern populations of Group 3 
all have very low heterozygosity and all are virtually identical genetically. 
The northern populations of groups 6, 12, and 13 have only a slightly 
lower than average variability. Thus only Group 3 shows the type of 
fixation of variable loci that might have occurred in the source popula- 
tions at the northern margin of the range during the glacial maxima. 
Some of these populations might have become very small, living as they 
did at the northern limit of survivable climate and habitat where genetic 
drift might have resulted in the fixation of variable loci. 

Because only one locus (Pt-3) is monomorphic among all samples, 95 
percent of the loci studied are genetically variable. 


Size Variation among Groups 


Duncan and Highton (1979) mention some of the problems associated 
with using adult size as a taxonomic character in salamanders. In spite 
of the expected difficulties, they found that there is considerable 
uniformity in maximum size within genetically differentiated popula- 
tions of the four species of the Plethodon glutinosus group in the Ouachita 
Mountains of Arkansas and Oklahoma. In the P. glutinosus complex, 
however, there appears to be an environmentally related pattern of 
geographic variation in maximum size. With the exception of P. aureolus 
and P. kentucki (small Appalachian species sympatric with larger 
members of the complex), size is correlated with physiographic 
province. 

Small-sized populations (maximum size usually under 74 mm snout 
to anterior angle of the vent) occur in the Coastal Plain physiographic 
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province (Table 6). In populations in the other upland physiographic 
provinces, the maximum size is usually larger. These results agree with 
the data in Highton (1962a). In the three groups that are distributed 
mostly in the coastal plain but also occur in other adjacent provinces 
(groups 3, 5, and 8), the maximum adult size of the non-coastal plain 
populations is greater than that of the coastal plain populations. This 
pattern might be the result of ecophenotypic variation, perhaps related 
to nutrition. If genetically based, however, selection might favor 
different adult sizes in the different soil types of the two areas. The 
influence of the sandy soils of the coastal plain versus the rocky soils of 
the other provinces in relation to maximum burrow size should be 
investigated. 


Table 6. Variation in the maximum size (snout to anterior angle of vent in life) of 
members of the P. glutinosus complex. 


Range of 
Group Species N maximum sizes Mean 
(mm) 
Coastal Plain populations 
3 P. chlorobryonis 9 59-73 68.1 
4 P. variolatus 7 63-72 66.9 
5 P. ocmulgee 4 64-70 65.8 
6 P. glutinosus 1 72 72.0 
8 P. mississippi 10 64-78 ait 
9 P. kisatchie 3 71-73 72.0 
11 P. grobmani 10 63-69 66.2 
13 P. albagula 1 70 70.0 
14 P. savannah 1 69 69.0 
Non-Coastal Plain populations 
1 P. teyahalee 6 72-90 81.0 
2 P. chattahoochee 4 72-79 75.0 
3 P. chlorobryonis 3 77-81 79.0 
5 P. ocmulgee 1 81 81.0 
6 P. glutinosus 39 68-89 78.7 
7 P. kiamichi 2 70-83 76.5 
8 P. mississippi 1 84 84.0 
10 P. sequoyah j 78 78.0 
2 P. cylindraceus 13 76-88 81.4 
13 P. albagula 12 68-82 74.5 
15 P. aureolus 3 57-67 63.7 
16 P. kentucki 4 64-75 68.5 
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Further analysis of geographic size variation will become available 
when my on-going study of the life histories of all eastern species of the 
genus Plethodon is completed. 


Taxonomy 


The mean D values of comparisons among all 16 groups are shown in 
Table 3. Of the 120 comparisons, only 14 are 0.20 or less, 36 are between 
0.21 and 0.29, and 70 are 0.30 or higher. Thus there is a considerable 
amount of genetic divergence among most of the groups. Avise (1974) 
and Ayala (1975) indicate that local populations of species commonly 
have Nei genetic distances less than 0.1. Avise and Aquadro (1982) 
showed that there is considerable variation in genetic distances between 
congeneric species among classes of vertebrates with birds having the 
lowest and amphibians the highest. 

Genetic distances between subspecies vary considerably because this 
taxonomic category is used to indicate so many different relationships. 
Avise (1974) states that “in many cases, biochemical systematists have 
not been able to distinguish subspecies which have been described by 
classical systematic criteria.” Many of these subspecies were probably 
based on minor geographic morphological or color pattern differentia- 
tion and are not the result of long periods of evolutionary divergence. 
Some subspecies recognized by taxonomists, however, may have evolved 
during long periods of isolation or geographic differentiation and in 
these the amount of genetic divergence might be greater, as well as 
concordant with the geographic patterns of morphological variation. 
An extreme example is the designation by Drosophila geneticists of 
highly differentiated populations of the D. willistoni group as subspecies. 
These subspecies have a mean D = 0.23 (Ayala, 1975: Table VI). 
Subspecies in this group are defined by Dobzhansky et al. (1977: 192) as 
“geographically separated, showing only a trace or no ethological 
isolation, but producing sterile male hybrids in at least one of the 
reciprocal crosses” in the laboratory. Most modern taxonomists do not 
hold this concept of subspecies and usually recognize subspecies as 
differentiated, but naturally freely interbreeding, geographical repre- 
sentatives of a species. The mean D (0.23) between populations consid- 
ered by Ayala (1975) to be “semispecies” (also partially reproductively 
isolated forms) is the same as between his “subspecies,” indicating a 
similar amount of genetic divergence. Because his “subspecies” and 
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“semispecies’ differ only in relative amount of reproductive isolation, as 
tested in the laboratory, both should be considered species taxonom- 
ically. Neither should be considered a subspecies because they have 
already evolved a significant degree of reproductive isolation (sterility in 
male hybrids). 

In much of the work being done on genetic variation in vertebrate 
species (see Smith et al., 1982, for a bibliography), evolutionary biolo- 
gists do not appear willing to make the obvious taxonomic conclusions 
that follow from their biochemical data. For example, Larson (1983) 
found that one population of Bolitoglossa occidentalis (sample 17) was 
genetically much closer to another species (B. rufescens) than it was to the 
three other samples of B. occidentalis. Moreover, his seven samples of 
B. rufescens were obviously divided into two genetically differentiated 
groups. Although Larson recognized that the present taxonomic 
arrangement of these populations into two species (based on the 
presence or absence of maxillary teeth) must be in error, he did not make 
the taxonomic allocations of populations indicated by his genetic data, in 
spite of the fact that his results were based on 29-30 genetic loci and the 
present taxonomy is based on only one morphological character that he 
states “is not a reliable indicator of species boundaries.” It should be 
pointed out, however, that the B. rufescens group may contain several 
species whose precise boundaries cannot be determined from the 
available data. Larson’s data indicate that his populations 17—21 have 
been separated from populations 22—27 for approximately 12 million 
years. Surely these two groups (and perhaps additional subgroups 
within both groups) are separate species and should not continue to be 
taxonomically arranged in the old erroneous fashion (populations 17, 
22, 26, and 27 as B. occidentalis and the rest as B. rufescens). 

Why are many biologists reluctant to make the taxonomic revisions 
that are indicated by the genetic data obtained from their biochemical 
research? The reasons are probably among the following: 

(1) Many biologists doing biochemical research on genetic variation 
are not systematists and do not understand or care to employ the 
procedures used by taxonomists to incorporate the conclusions of 
studies on the relationships of organisms into the formal system of 
nomenclature. Some journals do not accept species descriptions and 
these must therefore appear in separate papers from the biochemical 
results. 

(2) In the absence of direct information on reproductive isolation, 
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some workers are hesitant to recognize taxa that are differentiated 
genetically but not morphologically (Wake, 1981), or in which the 
biochemical data disagree with the morphological data, since the 
present taxonomic system is based largely on morphology. 

(3) Many nontaxonomists have too high a regard for the “current 
taxonomy.” Were they more familiar with the meager justification that 
sometimes supports present taxonomic allocations, they would be far 
less likely to accept uncritically the judgment of the last taxonomist who 
happened to revise the group with which they are concerned. 

(4) Genetic differentiation in structural genes (the type usually stud- 
ied by current biochemical techniques) is much more dependent on 
time since populations have shared a common ancestor than is the rate 
of change in morphology. Thus the most accurate phylogenetic recon- 
structions should be based on structural gene evolution. Surprisingly, 
some cladists (e.g., Farris et al., 1982) seem to be the most negative about 
the value of biochemical data in phylogenetic reconstruction and 
appear to be more interested in using Hennigian methods for mor- 
phological characters (well known to evolve at different rates in different 
evolutionary lines). They reject the very characters (products of struc- 
tural genes) that should provide the best data for estimating the real 
phylogeny of a group since they apparently evolve at much more 
constant rates than morphological characters. 

(5) Taxonomists have long recognized that the amount of morphologi- 
cal differentiation at the species level is less important in determining 
taxonomic status than is distinctness (Mayr, 1957). Thus taxonomists are 
conditioned to look for distinctive features separating sympatric popu- 
lations as indications of reproductive isolation rather than using the 
amount of morphological differentiation for that purpose. Systematists 
not familiar with the new biochemical techniques are understandably 
cautious concerning genetic distance data that cannot be used as 
diagnostic or key characters in the same way as standard taxonomic 
characters. However, electrophoretic data often provide information 
not only on amount of differentiation but also on whether or not there 
are fixed genetic differences between sympatric forms (indicating re- 
productive isolation in nature). 

(6) The taxonomic recognition of isolated allopatric populations has 
always been difficult for the taxonomist because the definitive criterion 
of reproductive isolation in nature cannot be applied. I have argued (zn 
Wake, 1981: 261—262) that in salamanders, which have so few variable 
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systematic characters and may evolve morphologically at extremely slow 
rates, genetic distance data may be more useful in determining the 
taxonomic status of allopatric populations than morphological data. 
Wake (1981) argues that chaos would result if the amount of genetic 
divergence, as measured by allozyme data, were used as a basis for 
making taxonomic decisions on the status of allopatric populations. He 
also feels the same way about the subdivision of continuously distributed 
populations into taxonomically recognized units on the basis of genetic 
distance data (personal communication). 

Patterns similar to those seen in the P. glutinosus complex have been 
found in western plethodontid salamanders of the genus Batrachoseps 
(Yanev, 1978, 1980). Yanev discovered several parapatric genetically 
differentiated forms that have variable amounts of gene exchange in 
contact zones where their ranges meet or overlap. She recognized these 
forms as semispecies although she (Yanev, 1980) applies to them the 
trinomen (nomenclaturally indicating that they are subspecies). How- 
ever, I agree that the semispecies seems to best indicate the taxonomic 
relationships of these populations and it is the arrangement I suggest 
for groups 1-14. Groups 15 and 16 (P. aureolus and P. kentucki) are widely 
sympatric with other groups and are not semispecies. Mayr (1963) 
defines semispecies as “the component species of a superspecies; also 
populations that have acquired some, but not yet all, attributes of 
species rank; borderline cases between species and subspecies.” 

The taxonomic recognition of parapatric or allopatric genetically 
differentiated forms is often equivocal for several reasons: 

(1) There is no established criterion of a specific amount of genetic 
differentiation that is always associated with the development of repro- 
ductive isolation. Some workers, for example Wake (1981), believe that 
there is as yet no convincing evidence that build-up of genetic distance, 
as measured by allozyme variation, should alone be used to determine 
species status, which he believes primarily should be based on mor- 
phological data. Although Wake correctly emphasizes the value of 
electrophoretic data in properly determining the real taxonomic rela- 
tionships in some instances (e.g., when two morphotypes of a popula- 
tion are genetically identical in allele frequences at all loci and are clearly 
interbreeding, or when two sibling sympatric forms have fixed genetic 
differences and are therefore not interbreeding), he (personal commu- 
nication) does object to making species level determinations of tax- 
onomic status solely on the basis of allozyme data, particularly when the 
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populations are allopatric. Baverstock et al. (1977), Thorpe (1982), and 
Bullini (1983), however, believe that estimates of genetic distance are 
probably the best criteria for making taxonomic decisions in doubtful 
cases. Borderline cases could still be very difficult because the confi- 
dence intervals of Nei distances tend to be quite large unless a great 
many more loci are used than in most current electrophoretic studies. 
In addition, Uzzell and Wake (personal communication) believe there is 
a substantial difference in comparisons in which there are merely 
frequency differences between two populations compared to cases 
where there are fixed differences, even though the D values in the two 
examples may be the same. The particular mix of fast- and slow- 
evolving loci (Sarich, 1977) used in a study also has a considerable effect 
on the magnitude of the D values obtained. 

(2) The treatment of isolated allopatric populations has always been a 
difficult problem for the taxonomist. Mayr (1942) indicates that the 
designation of such populations as species or subspecies is often subjec- 
tive because the extent of reproductive isolation cannot usually be 
determined. He suggests that a measure of the amount of morphologi- 
cal divergence that usually is associated with reproductive isolation 
between good species in each taxonomic group can be estimated by 
comparing closely related sympatric species of the same group. Brown 
and Wilson (1956) urge caution in using this measure for allopatric 
forms because of the likelihood of morphological character displace- 
ment in closely related sympatric species. I suggest that genetic 
distance data may be superior to morphological divergence in aiding 
the taxonomist to make better species level determinations in such 
situations. 

(3) In many of the 16 P. glutinosus groups there are no known 
diagnostic morphological characters that may be used to identify either 
living or preserved specimens. However, museum specimens of the 
other previously recognized species in the group (yonahlossee, jordani, 
ouachitae, caddoensis, and fourchensis) are also difficult to distinguish 
after their diagnostic color pigments have been lost in preservatives. 
Although the same is true of many of the 16 groups of the P. glutinosus 
complex, that should not preclude their taxonomic recognition. No 
careful multivariate morphometric study has yet been made on any of 
these species of large plethodons. The advent of electrophoretic 
methods of analysis now make it possible not only to identify species but 
to determine the geographic subgroup within each species with consid- 
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erable accuracy, but effective use of the method requires living individu- 
als of all of the forms or a complete frozen tissue collection. 

(4) The taxonomic recognition of forms that interbreed in narrow 
contact or overlap zones has long been a difficult one for taxonomists. 
Although evolutionarily these populations are intermediate between 
subspecies (freely interbreeding populations of the same species that 
differ morphologically) and species (reproductively isolated popula- 
tions), their taxonomic status is often difficult to decide. It is expected 
that all pairs of gradually diverging populations would be in an 
intermediate stage of differentiation at some time during their evolu- 
tion. If a number of populations of the ancestral P. glutinosus were all 
isolated at approximately the same time (as suggested above), then 
perhaps all might be in this intermediate period at the same time. The 
category semispecies has been suggested for such parapatric forms that 
have acquired some but not all the characteristics of species or are 
borderline cases between species and subspecies (see above). Amadon 
(1966) distinguishes between allospecies (parapatric species comprising 
a superspecies) and semispecies (populations that still exchange genes 
but not as freely as among conspecific populations and thus have only 
partially completed the process of speciation). Mayr (1963) states that 
“semispecies are a special kind of species, not a category different from 
a species.” Others, such as Hall and Selander (1973), Dobzhansky et al. 
(1977), and Yanev (1980), call such forms semispecies but refrain from 
naming them. Hall and Selander (1973) and Barton and Hewitt (1983) 
argue that strong selection against hybrids or backcross progeny in 
contact zones is essentially a barrier to gene flow between such para- 
patric forms, causing these zones to be very narrow. On this basis they 
regard the forms as species because they are essentially reproductively 
isolated. At the opposite extreme, Key (1981) argues that reproductive 
isolation should be defined as the relationship between two forms that 
do not ever hybridize in nature or whose F, hybrids leave no fertile 
progeny. 

Thus in this one complex, previously recognized as a single taxonomic 
species, we find many of the problems that have troubled taxonomists in 
defining the limits of species. In previous papers (Highton and 
MacGregor, 1983; Highton, 1984), four of the groups (1, 6, 15, and 16) 
have been recognized taxonomically as species because they are sym- 
patric with and reproductively isolated from other groups. The 
remaining 12 groups replace groups | and 6 and each other geograph- 
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ically. In some narrow contact zones there is evidence of gene exchange, 
whereas in others, such as those separated by large rivers, fixed differ- 
ences indicate no gene exchange for a considerable amount of time. No 
data are yet available for some contact zones, but the amount of genetic 
divergence is large. In still others there appears to be considerable gene 
exchange as indicated by several populations with intermediate allele 
frequencies. Genetic variation across contact zones is now being studied 
by the writer. 

There can be no question about the recognition of P. kentucki and 
P. aureolus as distinct species because they are reproductively isolated 
forms sympatric with other members of the complex. I feel that the 
remaining groups should be recognized as semispecies. Some taxono- 
mists would probably prefer to recognize them as subspecies and this 
option is always open to those who believe that the amount of genetic 
divergence is not sufficient to recognize these forms as species. For 
those who object to using allozyme data in making taxonomic judge- 
ments, groups I-14 can be referred to collectively as the “P. glutinosus 
complex.” 

My taxonomic treatment of the Plethodon glutinosus complex is the 
same as that adopted for leopard frogs of the Rana pipiens complex. 
The latter is also comprised of many closely related parapatric and 
allopatric forms that sometimes hybridize in narrow contact zones. 
The amount of genetic differentiation, as indicated by call differences 
and partial reproductive isolation when crossed in the laboratory, is 
sufficient to support their recognition as species rather than sub- 
species (Pace, 1974; Frost and Platz, 1983; Hillis et al., 1983; Platz 
and Frost, 1984; and references therein). Unpublished data on elec- 
trophoretic genetic comparisons among the various species of the R. 
pipiens complex by Sage (personal communication) indicate that their 
genetic distances are in the same range as those of the P. glutinosus 
complex. 

Whether or not all of the forms of the P. glutinosus complex now have 
been recognized remains to be determined. As mentioned above, 
P. kentucki has highly differentiated populations that are not analyzed 
in this paper, and it may be comprised of more than one form. 
Moreover, in spite of the extensive sampling, some types recognized 
here could have been missed because of their very small ranges (e.g., 
groups 7, 10, and 14). Further work could well discover additional 
forms. 
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Some of the following synonymies are incomplete in that the first 
reference to each species as Plethodon glutinosus is usually omitted, 
except in cases where salamanders from the original locality have been 
studied electrophoretically. If morphometric studies eventually pro- 
vide characters to distinguish between the species of the Plethodon 
glutinosus complex recognized in this paper, then it should be pos- 
sible to assign individuals mentioned in the literature to the proper 
form, provided investigators placed voucher specimens in museum 
collections. 


Group 1. 
Plethodon teyahalee Hairston 
Southern Appalachian Slimy Salamander 


Plethodon glutinosus (Green): Brimley (1912) (part). Highton (1970) 
(part). 

Plethodon jordan: teyahalee Hairston (1950:269). 

Plethodon jordan Blatchley: Highton (1962b). 

Plethodon (glutinosus) glutinosus (Green): Bishop (1941) (part). 
P(lethodon). teyahalee Hairston: Highton (1984). 


Holotype: UMMZ 100807, an adult male collected 23 August 1949 by 
Nelson G. Hairston, on Teyahalee Bald (= Johanna Bald), at an eleva- 
tion of 1,380 m in the Snowbird Mountains, Graham-Cherokee county 
line, North Carolina. 


Diagnosis: A large, light-chinned species with very small white dorsal 
spots, reduced lateral white spotting, and often with small red spots on 
the legs. The unique combination of genetic alleles that distinguishes 
P. teyahalee from other species of the P. glutinosus group is shown in Table 
4. Pgiallele c and Trf allele a are characteristic of P. teyahalee populations 
but are usually rare or absent in the other species. 


Distribution: West of the French Broad River in the Blue Ridge physio- 
graphic province of southwestern North Carolina and in immediately 
adjacent Tennessee. It also occurs in northern Rabun County, Georgia, 
and in Oconee, Pickens, Anderson, and Abbeville counties, South 
Carolina (Fig. 8). 
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Fig. 8. Distribution of P. teyahalee (Group 1) in the southern Appalachian 
Mountains based on samples identified to species electrophoretically. 


Group 2. 
Plethodon chattahoochee Highton, new species 
Chattahoochee Slimy Salamander 


Plethodon shermani Brimley: Howell (1909). 

Plethodon glutinosus (Green): Bailey (1937) (part). Highton (1970) (part). 
Plethodon shermani rabunensis Pope and Hairston: Hairston and Pope 
(1948) (part). 

Plethodon jordani rabunensis Pope and Hairston: Hairston (1950) (part). 
Plethodon jordani Blatchley: Highton (1962b) (part). 

Plethodon glutinosus glutinosus (Green): Highton (1962b) (part). 


Holotype: USNM 168527, an adult male collected 22 July 1961 by 
Richard Highton and Thomas Savage, at locality 9 (Fig. 1, Appendix 1), 
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0.3 km east of the top of Brasstown Bald, at an elevation of 1,353 m, 
‘Towns County, Georgia. 


Allotype: USNM 168518, an adult female, same data as holotype. 
Paratypes: USNM 168519-26, 168528-35; same data as holotype. 


Other material: Specimens from localities indicated in Figure 9 have been 
identified genetically as P. chattahoochee. 


Diagnosis: A large, light-chinned species with little or no dorsal spotting 
(except for populations in the extreme western part of its range) and 
abundant white or yellow lateral spotting. The unique combination of 
genetic alleles that distinguishes P. chattahoochee from other species of 
the P. glutinosus group is shown in Table 4. Allele c at the Pt-2 locus is 
unique to this species. 


Description of Holotype: After preservation: lengths from tip of snout to 
anterior angle of vent, 59 mm (in life, 62 mm); to posterior angle of vent, 
64 mm; total length, 115 mm. Head length (snout to gular fold) 14 mm; 
head width at widest point, 9 mm; length of front limbs, 14 mm, and of 
hind limbs, 15 mm. Vomerine teeth number 7 on both sides. 16 costal 
grooves (equivalent to 17 trunk vertebrae). Coloration in life: no white 
or brassy spotting on the black dorsum; numerous yellow spots on sides 
and a few small yellow spots on chin and underside of legs; few tiny red 
spots on dorsal surface of front limbs; venter lighter than ground color 
of back and sides; three small yellow spots on anterior venter; a few tiny 
brassy flecks on both eyelids and one small white spot on right eyelid; a 
trace of brassy flecking in iris of both eyes. Chin lighter than venter and 
with a prominent, round mental gland. Many tiny spots (hedonic 
glands) on ventral surface of body and tail. 


Description of Allotype: After preservation: length from tip of snout to 
anterior angle of vent, 60 mm (in life, 64 mm), to posterior angle of vent, 
65 mm; total length, 132 mm; head length, 14 mm; head width, 9 mm; 
length of front limbs, 13 mm, and of hind limbs, 16 mm. Vomerine teeth 
number 8 on the right and 9 on the left. 16 costal grooves. Coloration in 
life: similar to holotype except no red spots on legs, no yellow spots on 
venter, and no brassy flecks in iris. 
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Variation in the Type Series: In life, 9 of 18 specimens in the type series had 
small red spots on anterior legs; only the smallest (22 mm snout-anterior 
angle of vent) had a red spot on the dorsum. None had any red pigment 
on the hind limbs, but one had one tiny red spot on one cheek. The 
juvenile had several brassy flecks on the anterior dorsum, but only two 
of the larger animals had traces of dorsal brassy flecking (one only on 
the head, the other only on the dorsum). All had yellow spotting on the 
sides and undersides of the limbs. Twelve of the 18 had yellow spotting 
on the chin and 8 also had similar spotting on the anterior venter. All 
but the allotype had brassy flecking in the irises and all but one had 
brassy flecking on the eyelids. 


Distribution: Most of the Blue Ridge physiographic province of northern 
Georgia (Fig. 9). It also occurs in southeastern Cherokee County, North 
Carolina. The distribution of this species approximately coincides with 
the limits of the Chattahoochee National Forest, after which the species 
is named. 


Group 3. 
Plethodon chlorobryonis Mittleman, new status 
Atlantic Coastal Slimy Salamander 


Plethodon glutinosus chlorobryonis Mittleman (1951:108). 


Holotype: USNM 129933, an adult male collected 10 November 1950 by 
Myron B. Mittleman and C. B. Goodstein, at or near locality 14 (Fig. 1, 
Appendix 1), 13 miles north of New Bern, Craven County, North Carolina. 


Diagnosis: A small-sized species (except that non-coastal plain popula- 
tions tend to be larger) with very small, slightly brassy, white dorsal spots 
and abundant lateral white or yellow spotting. The unique combination 
of genetic alleles that distinguishes P. chlorobryonis from other species of 
the P. glutinosus group is shown in Table 4. 


Distribution: Coastal Plain physiographic province of southeastern 
Virginia, North Carolina, and northeastern South Carolina. It enters 
the Piedmont physiographic province in southeastern Virginia, central 
and western South Carolina, and the Blue Ridge physiographic prov- 
ince in northeastern Georgia (Fig. 10). 
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Fig. 10. Distribution of of P. chlorobryonis (Group 3) in the middle Atlantic states 
based on samples identified to species electrophoretically. 
Group 4. 


Plethodon variolatus (Gilliams), new combination 


South Carolina Coastal Slimy Salamander 


Salamandra variolata Gilliams (1818:460). 
Plethodon variolosum (Gilliams). Duméril, Bibron, and Duméril (1854). 


Holotype: Not known to exist. The type locality is the “southern states.” 
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Schmidt (1953) restricted the type locality to the vicinity of Charleston, 
South Carolina, but Neill (1957) criticized Schmidt’s restriction of this 
and other type localities. To remove future question as to which species 
this name should apply, I designate UNSM 267104 as the neotype of 
Salamandra variolata, an adult male collected 27 March 1986 by David E. 
Carr, from locality 27 (Fig. 1, Appendix 1), at an elevation of 6 m, 
Beechtree Recreation Area, Berkeley County, South Carolina. 


Description of Neotype: In life: length from tip of snout to anterior angle 
of vent, 57 mm, and to posterior angle of vent, 61 mm; total length, 121 
mm; head length, 13 mm; head width, 10 mm; length of front limbs, 13 
mm, and of hind limbs, 14 mm. (After preservation, snout-vent lengths 
52 and 56 mm and total length 114 mm.) Vomerine teeth number 8 on 
the right and 7 on the left. 16 costal grooves. A few small white dorsal 
spots, each with scattered brassy flecking, especially at edges of spots. 
Abundant large yellow spots on sides; few scattered yellow spots on chin 
and venter. Venter lighter than black ground color of dorsum; chin 
much lighter than venter and with a small round mental gland. 


Diagnosis: A small species with small to medium-sized dorsal brassy 
white spots and abundant lateral white or yellow spotting. Populations 
from Jasper County, South Carolina, lack the lateral and dorsal spot- 
ting, as do occasional specimens from other populations. The unique 
combination of genetic alleles that distinguishes P. variolatus from other 
species of the P. glutinosus group is shown in Table 4. 


Distribution: Southern Coastal Plain physiographic province of South 
Carolina and extreme southeastern Georgia. 


Group 5. 
Plethodon ocmulgee Highton, new species 
Central Georgia Slimy Salamander 


Holotype: USNM 257426, an adult male collected 8 November 1976 by S. 
Blair Hedges and Richard Highton, at locality 32 (Fig. 1, Appendix 1), 
Little Ocmulgee State Park, at an elevation of 49 m, Wheeler County, 
Georgia. 


Allotype: USNM 257427, an adult female, same data as holotype. 
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Paratypes: USNM 257428-257464, same data as holotype. 


Diagnosis: A small-sized species (except for Piedmont populations) with 
few small brassy dorsal spots and a moderate amount of lateral white 
spotting. The unique combination of genetic alleles that distinguishes 
P. ocmulgee from other species of the P. glutinosus group is shown in 
Table 4. Allele a of the Gdh locus, rare or absent from all other species, 
is always common in P. ocmulgee. 


Description of Holotype: After preservation: length from tip of snout to 
anterior angle of vent, 52 mm (in life, 56 mm) and to posterior angle of 
vent, 57 mm; total length, 103 mm. Head length (snout to gular fold), 
13 mm; head width at widest point, 8 mm; length of front limbs, 12 mm, 
and of hind limbs, 13 mm. Vomerine teeth number 4 on the right and 7 
on the left. 16 costal grooves. Few small brassy and white dorsal spots in 
life and numerous white spots on sides; venter almost as dark as black 
ground color of back and sides, chin slightly lighter than venter; many 
tiny spots (hedonic glands) on ventral surface of the body and tail; a 
small round mental gland. 


Description of Allotype: After preservation: length from tip of snout to 
anterior angle of vent, 59 mm (in life, 62 mm) and to posterior angle of 
vent, 64 mm; total length, 119 mm; head length, 15 mm; head width, 9 
mm; length of front limbs, 13 mm on the right and 14 mm on the left; 
and of hind limbs, 15 mm. Vomerine teeth number 8 on each side. 16 
costal grooves. Coloration in life similar to holotype. 


Distribution: The upper Coastal Plain and adjacent Piedmont physio- 
graphic provinces of central Georgia. Much of its range is in the 
Ocmulgee River drainage, after which the species is named. 


Group 6. 
Plethodon glutinosus (Green) 
Northern Slimy Salamander 


Salamandra glutinosa Green (1818:357). 
Plethodon glutinosus (Green): Tschudi (1838). 
Plethodon glutinosum (Green): Gray (1850). 
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Cylindrosoma glutinosum (Green): Dumeéril, Bibron, and Duméril (1854). 
Salamandra melanoleuca Wied (1865:130—131). Type locality: Nazareth, 
Pennsylvania. 

Plethodon glutinosus glutinosus (Green): Dunn (1920). Bishop (1943). 
Grobman (1944). 


Holotype: Dunn (1926) states that the holotype is not known to exist, but 
the type locality is probably Princeton, New Jersey. 


Diagnosis: A large species with large brassy colored dorsal spots and 
moderately abundant lateral white or yellow spotting. It differs from 
sympatric P. kentucki in having a darker chin, smaller mental gland in 
adult males, and larger, more brassy dorsal spots. It is not detectably 
different morphologically from groups 7, 10, and 13. The unique 
combination of genetic alleles that distinguishes P. glutinosus from other 
species of the P. glutinosus group is shown in Table 4. 


Distribution: Southwestern Connecticut west to southern Illinois and 
south through West Virginia, western Virginia, Kentucky and Tennes- 
see to eastern Alabama and northwestern Georgia (Fig. 11). 


Group 7. 
Plethodon kiamichi Highton, new species 
Kiamichi Slimy Salamander 


Holotype: USNM 257314, an adult male collected 5 June 1973 by Wayne 
Garber, Richard Highton, and James Hook, at locality 76 (Fig. 1, 
Appendix 1), Round Mountain, at an elevation of 640 m, LeFlore 
County, Oklahoma. 


Allotype: USNM 257315, an adult female, same data as the holotype. 
Paratypes: USNM 257316-257347, same data as the holotype. 


Other material: Specimens from localities indicated in Figure 12 have 
been identified genetically as P. kiamichi. 


Diagnosis: A large species with large brassy dorsal spots and moderately 
abundant lateral white or yellow spotting. It is not detectably different 
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morphologically from groups 6, 10, and 13. The unique combination of 
genetic alleles that distinguishes P. kiamichi from other species of the P. 
glutinosus group is shown in Table 4. Got-1 allele 6, common in P 
kiamichi, is absent from the other species of the P. glutinosus complex, 
but is also common in P. ouachitae (Table 4). 


“‘W1HO 


Latimer 


Kiamichi Mountain 


Pushmataha 


McCurtain 
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Fig. 12. Distribution of P. kiamichi (Group 7) on Kiamichi Mountain, Arkansas 
and Oklahoma, based on samples identified to species electrophoretically. 
Counties are indicated. 


Description of Holotype: After preservation: length from tip of snout to 
anterior angle of vent, 65 mm (in life, 69 mm), and to posterior angle of 
vent, 70 mm; total length, 141 mm; head length (snout to gular fold), 15 
mm; head width at widest point, 10 mm; length of front limbs, 15 mm, 
and of hind limbs, 16 (right) and 15 mm (left). Vomerine teeth number 
13 on the right and 11 on the left. 16 costal grooves. Many small white 
dorsal spots with associated brassy flecking in life; numerous white 
spots on sides, chin, and legs; venter dark as ground color of back and 
sides and with few white spots; chin as dark as venter with a prominent, 
large, oval-shaped mental gland. Many tiny spots (hedonic glands) on 
ventral surface of body and tail. Few tiny brassy flecks and white spots 
on the eyelids in life. 


Description of Allotype: After preservation: length from tip of snout to 
anterior angle of vent, 79 mm (in life, 83 mm), and to posterior angle of 
vent, 85 mm; total length, 160 mm; head length, 17 mm; head width, 11 
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mm; length of front limbs, 15 mm, and of hind limbs, 17 mm. Vomerine 
teeth number 9 on the right and 8 on the left. 17 costal grooves. 
Coloration in life similar to holotype except for absence of ventral 
hedonic glands. 


Distribution: Known only from Round and Kiamichi mountains in Polk 
County, Arkansas, and LeFlore County, Oklahoma (Fig. 12). This 
species is named for Kiamichi Mountain. 


Group 8. 
Plethodon mississippi Highton, new species 
Mississippi Slimy Salamander 


Holotype: USNM 257388, an adult male collected 18 January 1986 by 
David E. Carr and Richard Highton, at locality 79 (Fig. 1, Appendix 1), 
Tishomingo State Park, at an elevation of 177 m, Tishomingo County, 
Mississippl. 


Allotype: USNM 257389, an adult female, same data as the holotype. 
Paratypes: USNM 257390-257425, same locality as the holotype. 


Other material: Specimens from localities indicated in Figure 13 have 
been identified genetically as P. massissippr. 


Diagnosis: A small species (except for populations in the Appalachian 
Plateau of northern Alabama that reach a larger size) with large brassy 
dorsal spots and abundant lateral white or yellow spotting. Mor- 
phologically, it is not detectably different from Group 11. The unique 
combination of genetic alleles that distinguishes P. mississeppi from other 
species of the P. glutinosus group is shown in Table 4. 


Description of Holotype: In life: length from tip of snout to anterior angle 
of vent, 65 mm, and to posterior angle of vent, 69 mm; total length, 
137 mm; head length, 15 mm; head width, 10 mm; length of front 
limbs, 15 mm, and of hind limbs, 16 mm. Vomerine teeth number 10 on 
the right and 11 on the left. 16 costal grooves. Many small white spots 
with slight brassy flecking on black dorsum; numerous yellow spots on 
sides; few small white spots on chin and legs; belly as dark as ground 
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color of back and sides and chin slightly lighter than belly with a round 
mental gland. Many tiny spots (hedonic glands) on ventral surface of 
body and tail. Few tiny brassy flecks and small white spots on eyelids. 
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Fig. 13. Distribution of P. massissippi (Group 8) in the south-central United States 
based on samples identified to species electrophoretically. 


Description of Allotype: After preservation: length from tip of snout to 
anterior angle of vent, 66 mm (in life, 71 mm), and to posterior angle of 
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vent, 71 mm; total length, 142 mm; head length, 15 mm; head width, 10 
mm; length of front limbs, 14 mm; and of hind limbs, 16 mm. Vomerine 
teeth number 7 on the right and 9 on the left. 16 costal grooves. 


Distribution: Western Alabama, Mississippi, Florida Parishes of south- 
eastern Louisiana, western Tennessee, and western Kentucky (Fig. 13). 
The species is named for the State of Mississippi. 


Group 9. 
Plethodon kisatchie Highton, new species 
Louisiana Slimy Salamander 


Holotype: USNM 257348, an adult male collected 27 January 1971 by 
Richard Highton, James Hook, George Kramer, Mark Mello, Donald C. 
Morizot, and David Walter, at locality 90 (Fig. 1, Appendix 1), along 
Indian Creek, at an elevation of 30 m, Grant Parish, Louisiana. 


Allotype: USNM 257349, an adult female, from locality 89 (Fig. 1, 
Appendix 1), Winn Parish, Louisiana, collected on 16 December 1973 
by Richard Highton and Scott Highton. 


Paratypes: USNM 257350-257355, same locality as the holotype; USNM 
257356- 257387, same locality as the allotype. 


Diagnosis: A small species with very large brassy dorsal spots. The 
unique combination of genetic alleles that distinguishes P. kisatchie from 
other species of the P. glutinosus group is shown in Table 4. Est allele / is 
usually rare or absent in the other species. 


Description of Holotype: After preservation: length from tip of snout to 
anterior angle of vent, 52 mm (in life, 56 mm), and to posterior angle of 
vent, 56 mm; total length, 115 mm; head length, 13 mm; head width, 8 
mm; length of front limbs, 12 mm, and of hind limbs, 13 mm. Vomerine 
teeth number 7 on the right and 8 on the left. 16 costal grooves. 
Coloration in life: many large brassy white spots on black dorsum; 
numerous yellow spots on sides and legs but none on chin and venter; 
venter slightly lighter than ground color of back and sides; chin slightly 
lighter than venter; chin with small round mental gland; many tiny 
spots (hedonic glands) on ventral surface of body and tail; a few tiny 
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brassy flecks and small white spots on eyelids. Areas free of black 
pigmentation along anterior margin of gular fold, ventral margin of 
limb insertions, and front limb joints. 


Description of Allotype: After preservation: length from tip of snout to 
anterior angle of vent, 58 mm (in life, 60 mm), and to posterior angle of 
vent, 63 mm; total length, 125 mm; head length, 14 mm; head width, 9 
mm; length of front limbs, 12 mm, and of hind limbs, 14 mm. Vomerine 
teeth number 6 on the right and 5 on the left. 16 costal grooves. 
Coloration in life: similar to holotype except no hedonic glands on 
venter and reduced pigmentation (not absence) at hind limb insertions; 
a few small yellow spots on the chin and belly. 


Distribution: From central Louisiana north to southern Arkansas. The 
species is named for the Kisatchie National Forest. 


Group 10. 
Plethodon sequoyah Highton, new species 
Southeastern Oklahoma Slimy Salamander 


Holotype: USNM 257485, an adult male collected 11 October 1976 by S. 
Blair Hedges and Richard Highton, at locality 91 (Fig. 1, Appendix 1), 
Beavers Bend State Park, at an elevation of 140 m, McCurtain County, 
Oklahoma. 


Allotype: USNM 257486, an adult female, same locality as the holotype, 
collected on 25 May 1978 by Richard Highton, Scott Highton and 
Jeffrey Streicher. 


Paratypes: USNM 257487-257521, same data as the allotype. 


Diagnosis: A large species with large brassy dorsal spots and moderately 
abundant white or yellow spotting. It is not detectably different mor- 
phologically from groups 6, 7, and 13. The unique combination of 
genetic alleles that distinguishes P. sequoyah from other species of the 
P. glutinosus group is shown in Table 4. Mdh-2 allele a is unique to 
P. sequoyah. 


Description of Holotype: After preservation: length from tip of snout to 
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anterior angle of vent, 58 mm (in life, 61 mm), and to posterior angle of 
vent, 63 mm; total length 133 mm; head length, 14 mm; head width, 9 
mm; length of front limbs, 14 mm; and of hind limbs, 15 mm. Vomerine 
teeth number 8 on the right and 9 on the left. 15 costal grooves. 
Numerous small white brassy dorsal spots; larger yellow lateral spots; 
chin with white spots; venter almost as dark as black ground color of 
back and sides, chin slightly lighter than venter with small, rounded, 
mental gland. 


Description of Allotype: After preservation: length from tip of snout to 
anterior angle of vent, 64 mm (in life, 69 mm), and to posterior angle of 
vent, 69 mm; total length, 136 mm; head length, 15 mm; head width, 10 
mm; length of front limbs, 14 mm, and of hind limbs, 15 mm. Vomerine 
teeth number 10 on the right and 9 on the left. 16 costal grooves. 
Coloration in life similar to holotype. 


Distribution: Known only from the type locality, Beavers Bend State 
Park, McCurtain County, Oklahoma. Named for the Oklahoma Indian 
Sequoyah, who devised the Cherokee alphabet. 


Group 11 
Plethodon grobmani Allen and Neill, new status 
Southeastern Slimy Salamander 


Plethodon glutinosus grobmani Allen and Neill (1949:112). 


Holotype: ERA-WTN 19220, an adult female collected 1 October 1949 by E. 
Ross Allen, Bobby Allen, and Wilfred T. Neill, at Half-mile Creek Swamp, 
about % mile northeast of Silver Springs, Marion County, Florida. 


Diagnosis: A small species with large brassy dorsal spots and abundant 
lateral white or yellow spotting. Morphologically, it is not detectably 
different from Group 8. The unique combination of genetic alleles that 
distinguishes P. grobmani from other species of the P. glutinosus group is 
shown in Table 4. 


Distribution: Southern Alabama and southern Georgia south to central 
Florida (Fig. 14). Topotypes are genetically similar to other populations 
from the Florida peninsula. 
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Fig. 14. Distribution of P. grobmani (Group 11) in the southeastern United States 
based on samples identified to species electrophoretically. 


Group 12 
Plethodon cylindraceus (Harlan), new combination 
White-Spotted Slimy Salamander 


Salamandra cylindracea Harlan (1825:156-157). 


Holotype: Not known to exist. The type locality is South Carolina, 
probably in the vicinity of Camden, although Schmidt (1953) restricted 
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the type locality to the vicinity of Charleston, South Carolina. Highton 
(1962b) pointed out that Camden is more likely the type locality than 
Charleston. Because Camden is in or near the contact zone between 
groups 3 and 12, I believe a neotype should be designated to avoid 
future confusion as to which name applies to this species. I therefore 
designate USNM 257522 as the neotype of Salamandra cylindracea, an 
adult female collected 31 March 1971 by Richard Highton and Donald 
C. Morizot, from locality 112 (Fig. 1, Appendix 1), at an elevation of 137 
m, Chester County, South Carolina. 


Description of Neotype: After preservation: length from tip of snout to 
anterior angle of vent, 70 mm (in life, 76 mm), and to posterior angle of 
vent, 74 mm; total length 147 mm; head length, 16 mm; head width, 11 
mm; length of front limbs, 14 mm; and of hind limbs, 18 mm. Vomerine 
teeth number 9 on the right and 7 on the left. 16 costal grooves. Dorsum 
and legs with small white spots in life and sides with larger yellow spots; 
a few small yellow spots on cheeks, chin, and venter. Venter and chin 
slightly lighter than black ground color of dorsum. 


Diagnosis: A large, light-chinned species with large white dorsal spots and 
moderately abundant lateral white spots. The unique combination of 
genetic alleles that distinguishes P. cylindraceus from other species of the 
P. glutinosus group is shown in Table 4. Mdh-] allele a, common in the north- 
ern populations of P. cylindraceus, is rare or absent in the other species. 


Distribution: The Piedmont and Blue Ridge physiographic provinces of 
Virginia and North Carolina west to the French Broad River and south 
to the northern Piedmont of South Carolina. It also occurs in parts of 
the Valley and Ridge physiographic province in western Virginia and 
extreme eastern West Virginia and in a small area of the Coastal Plain 
physiographic province of eastern Virginia (Fig. 15). 


Group 13. 
Plethodon albagula Grobman, new status 
Western Slimy Salamander 


Plethodon glutinosus albagula Grobman (1944:283). 


Holotype: CM 9652, an adult male collected 24 February 1935 by Wesley 
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Clanton, 20 miles north of San Antonio, Bexar County, Texas. 


Diagnosis: A large species with large brassy dorsal spots and moderately 
abundant lateral white or yellow spotting. Morphologically, it is not 
detectably different from groups 6, 7, and 10. The unique combination 
of genetic alleles that distinguishes P. albagula from other species of the 
P. glutinosus group is shown in Table 4. The null allele (h) in the Pt-2 


locus is found only in P. albagula. 
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Fig. 15. Distribution of P. cylindraceus (Group 12) in the middle Atlantic states 
based on samples identified to species electrophoretically. 


Distribution: Southern Missouri, the highlands of northern and western 
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Arkansas and eastern Oklahoma exclusive of the ranges of groups 7 and 
10; the Balcones Escarpment area of south-central Texas. Additional 
samples assigned genetically to this species are from Independence, 
Johnson, Logan, Polk, Scott, and Van Buren counties, Arkansas, and 
Warren County, Missouri (north of the Missouri River). 


Group 14. 
Plethodon savannah Highton, new species 
Eastern Georgia Slimy Salamander 


Holotype: USNM 257465, an adult male collected 17 March 1985 by 
Richard Highton and Hansjiirg Hotz, at locality 128 (Fig. 1, Appendix 
1), at an elevation of 101 m, Richmond County, Georgia. 


Allotype: USNM 257466, an adult female, same data as the holotype. 


Paratypes: USNM 257467-257480, same data as the holotype; USNM 
257481-257484 topotypes. 


Other material: Specimens from localities indicated in Figure 16 have 
been identified genetically as P. savannah. 


Diagnosis: A small species with very little brassy pigment in the 
white dorsal spots and with abundant lateral white or yellow spots. 
The unique combination of genetic alleles that distinguishes P. 
savannah from the other species of the P. glutinosus group is shown in 
Table 4. Idh-1 allele d is usually rare or absent in the other 
species. 


Description of Holotype: In life: length from tip of snout to anterior angle 
of vent, 67 mm, and to posterior angle of vent, 72 mm; total length, 139 
mm; head length, 15 mm; head width, 11 mm; length of front limbs, 13 
mm; and of hind limbs, 16 mm. Vomerine teeth number 9 on the right 
and 10 on the left. 16 costal grooves. Coloration in life: many white 
iridophore dorsal spots scattered on the black ground color; abundant 
white iridophore spotting on sides, dorsal surfaces of limbs, and moder- 
ate amount on cheeks, chin and sides of head; venter with few white 
spots. 
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Fig. 16. Distribution of P. savannah (Group 14) in eastern Georgia based on 
samples identified to species electrophoretically. Counties are indicated. 


Description of Allotype: In life: length from tip of snout to anterior angle 
of vent, 65 mm, and to posterior angle of vent, 70 mm; total length, 130 
mm; head length, 14 mm; head width, 10 mm; length of front limbs, 13 
mm; and of hind limbs, 16 mm. Vomerine teeth number 6 on the right 
and 7 on the left. 16 costal grooves. Coloration similar to holotype. 


Distribution: Known only from Burke, Jefferson, and Richmond coun- 
ties, Georgia (Fig. 16). The apparent eastern limit of the range is the 
Savannah River, after which the species is named. 
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Group 15. 
Plethodon aureolus Highton 
Tellico Salamander 


Plethodon aureolus Highton (1984:2). 


Holotype: USNM 238341, an adult male collected 30 June 1979 by 
Richard Highton and Jeffrey K. Streicher, at Farr Gap (locality 130; Fig. 
1; Appendix 1), at an elevation of 872 m, Monroe County, Tennessee. 


Diagnosis: The smallest species of the P. glutinosus complex with a light 
chin, abundant dorsal brassy spots, and abundant lateral white or yellow 
spotting. The unique combination of genetic alleles that distinguishes 
P. aureolus from other species of the P. glutinosus group is shown in Table 
4. Allele d at the Pt-2 locus is unique to P. aureolus. Allele fat the Ldh-m 
locus and Idh-2 allele d are usually rare or absent in the other species. 


Distribution: Between the Little Tennessee and Hiwassee rivers on the 
western slopes of the Unicoi Mountains in northeastern Polk and 
Monroe counties, Tennessee, and in western Graham and northwestern 
Cherokee counties, North Carolina (Highton, 1984: Fig. 1). 


Group 16. 
Plethodon kentucki Mittleman 
Cumberland Plateau Woodland Salamander 


Plethodon kentucki Mittleman (1951:105). Highton and MacGregor 
(1983). 
Plethodon jordani kentucki: Schmidt (1953). 


Holotype: USNM 129937, an adult male collected in August 1933 by W. 
Cornett, on Pine Mountain, at an elevation of about 610 m, Harlan 
County, Kentucky. 


Diagnosis: Geographically variable in size but generally smaller than 
sympatric P. glutinosus (Group 6). It usually has a lighter chin than 
sympatric P. glutinosus and its dorsal spots are smaller and possess less 
brassy pigment. The mental gland of adult males is larger than that of 
sympatric P. glutinosus. The unique combination of genetic alleles that 
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distinguishes P. kentucki from other species of the P. glutinosus group is 
shown in Table 4. Allele g of the Pt-2 locus is found only in P. kentucki, 
and Alb alleles a and 6, Trf allele m, Got-1 alleles a and d, and Got-2 allele 
d are usually rare or absent in the other species of the P. glutinosus 


complex. 


Distribution: The Cumberland Plateau of western West Virginia west of 
the New River, eastern Kentucky, southwestern Virginia, and north- 
eastern Tennessee (Fig. 17). 
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Fig. 17. Distribution of P. kentucki (Group 16) in the Cumberland Plateau based 
on samples identified to species electrophoretically. 


SUMMARY 


Genetic variation in 22 independent loci of Plethodon glutinosus was 
analyzed electrophoretically in 135 samples from populations collected 
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throughout its range. It was found that P. glutinosus is a complex 
comprised of at least 16 genetically divergent species and semispecies, 
most of which probably began differentiating during the late Miocene 
or Pliocene. One form, P. teyahalee, is probably of hybrid origin. 

The most different form, P. kentucki Mittleman, occurs in sympatry 
throughout the Cumberland Plateau with the northern, large, brassy- 
spotted species here recognized as P. glutinosus (Green). Two other 
highly differentiated species, P. teyahalee Hairston and P. aureolus 
Highton, are sympatric throughout the range of P. aureolus in south- 
eastern Tennessee and adjacent North Carolina. Both P. teyahalee and 
P. aureolus are known to be sympatric with P. glutinosus at a single 
locality in Polk County, Tennessee, without evidence of hybridization 
(Highton, 1984), but their ranges are otherwise parapatric to that of P. 
glutinosus. All forms except P. kentucki and P. aureolus replace each other 
geographically and are considered semispecies. Names are available for 
some of the species, but new names are provided for seven undescribed 
taxa. The problems associated with recognizing forms taxonomically 
that are primarily distinguished on the basis of allozyme data are 
discussed. 

The northernmost populations of P. chlorobryonis, P. glutinosus, and 
P. cylindraceus have much less geographic genetic variation than the 
more southern populations of most of the forms (an exception is the 
genetic similarity of the southern populations of P. massissippi). The 
genetic similarity of northern populations is probably due to the recent 
northward expansion of their ranges after the withdrawal of the Wis- 
consin glaciation during the last 15,000 years. 

Mean heterozygosity in populations varies from a high of 0.23 in two 
populations of P. chattahoochee, the most variable species, to a low of 0 in 
P. savannah and one population each of P. chlorobryonis and P. kisatchie. 
The overall mean heterozygosity of all 135 populations of the P. glutino- 
sus complex is 0.06. Of the loci evaluated electrophoretically, 21 of 22 
(95%) have detectable genetic variation. Differing rates of evolution are 
apparent at the loci evaluated. At some loci it is possible to suggest the 
ancestral genotype. 

The maximum size of coastal plain salamanders of all species is 
usually below 74 mm and that of noncoastal plain populations is usually 
above 74 mm. Exceptions are the small-sized species P. aureolus and P. 
kentucki. 
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Addendum 


Since this paper was prepared, it was discovered that a mistake in coding 
one locus in population 116 resulted in an error in its placement on the 
phenogram in Figure 1. When corrected, this, the only population of 
light-chinned Texas P. albagula, has an average Nei genetic distance to 
other populations of that species of 0.16, making it the most divergent 
population of P. albagula. Its relationship to other Texas populations is 
obviously in need of further study (see p. 19). 


Part II 
Immunological Analysis of Geographic 
Variation in Plasma Albumins 


George C. Maha, Linda R. Maxson, 
and Richard Highton 
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INTRODUCTION 


The quantitative immunological technique of micro-complement fixa- 
tion (MC’F) has been used extensively in studies of phylogenetic 
relationships among diverse amphibians (for example, Maxson et al., 
1979; Maxson, 1984; Maxson and Roberts, 1985; Hutchinson and 
Maxson, 1987). Micro-complement fixation studies of serum albumin 
(a 580 amino acid protein) demonstrate that MC’F offers a sufficiently 
rapid and accurate means of estimating sequence differentiation (Wil- 
son et al., 1977; Maxson and Maxson, 1986) to allow good phylogenetic 
inference over a long range of time (approximately 100 million years). 
Protein sequence information has a known genetic foundation and 
consequently can be analyzed phylogenetically and independently of 
any hypotheses about evolutionary rates. When evolutionary rates are 
unequal (and they generally do exhibit some variance), the resultant 
phylogenetic tree will display this fact of nonequal rates without any bias. 
Phylogenetic trees derived from molecular data do not group species 
together just because they are similar morphologically. Serum albumin, 
which evolves fairly rapidly, provides a data base such that the number 
of events between speciation nodes is large. Thus, each internodal 
lineage in a tree typically is defined by numerous events (generally ten 
or more). The number of detected changes along each of the lineages is 
a direct measure of the confidence in inferring that lineage. 
Micro-complement fixation analysis of proteins of known sequence 
has provided evidence that the immunological distance (ID) measured 
is a good estimate of the amino acid sequence difference between the 
two proteins compared (Champion et al., 1975; Wilson et al., 1977; 
Prager et al., 1978; Benjamin et al., 1984; Maxson and Maxson, 1986). 
For albumin, it is estimated that each ten ID units represents roughly 
ten amino acid substitutions in the albumins being compared (Maxson 
and Wilson, 1974; Maxson and Maxson, 1986). Using MC’F analyses of 
albumin evolution, single amino acid substitutions can be detected and 
comparisons between homologous albumins can be made until they 
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differ in their total sequence by 35-40% (Maxson and Maxson, 1979, 
1986). 

Because molecular evolution proceeds independently of morphologi- 
cal evolution (Wilson, 1985; Beverley and Wilson, 1982; Maxson et al., 
1979; Wilson et al., 1974), MC’F studies have uncovered cryptic species 
(Maxson et al., 1977; Maxson, 1981; Maha et al., 1983) and detected cases 
of convergent morphological evolution (Maxson and Wilson, 1974; 
Maxson, 1977; Maxson et al., 1982). Recent work has shown MC’F 
capable of identifying interspecific hybrids (Maxson et al., 1987). Thus, 
the application of MC’F as a sensitive probe of protein sequence evolu- 
tion has addressed many interesting questions in amphibian phylogeny 
and biogeography (for example Maxson and Roberts, 1985; Maxson et 
al., 1984; Maxson and Wake, 1981). 

Most of these earlier studies were concerned primarily with inter- 
specific rather than interpopulation comparisons. Several workers have 
shown that gel electrophoresis of many proteins, not MC’E, is the 
molecular technique of greatest sensitivity for studies of intraspecific 
protein variation (Bush and Kitto, 1978; Maxson and Maxson, 1979; 
Prager and Wilson, 1980). However, MC’F has been shown to be a 
sensitive detector of single amino acid differences between orthologous 
proteins (Cocks and Wilson, 1969; Maxson and Maxson, 1986). Indeed, 
some MC’F studies have compared allopatric populations of the same 
species and demonstrated some degree of albumin polymorphism 
between these populations (Maxson and Wilson, 1974; Maxson, 1981; 
Maha et al., 1983). 

While studying the phylogenetic relationships among the then recog- 
nized 26 species of Plethodon (Maxson et al., 1979), we prepared an 
antibody to albumin from a population of P. glutinosus from Ulster 
County, New York. Highton (1970, 1972) had shown that there were 
several locally differentiated populations of the P. glutinosus complex in 
the southern Appalachian Mountains and in the middle Atlantic states. 
Accordingly, we tested individuals from several populations that High- 
ton had designated as “brassy-spotted,” “white-spotted” or “Coastal 
Plain” types. Our preliminary survey demonstrated that MC’F compari- 
sons could distinguish P. glutinosus albumin alleles. The realization that 
we could readily measure albumin differentiation among these groups 
prompted us to begin a more extensive, simultaneous immunological 
and electrophoretic study of the P. glutinosus complex from throughout 
its range. 
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MATERIALS AND METHODS 


Antisera to purified albumin of P. glutinosus from Ulster County, New 
York, was available from the earlier study of Plethodon (Maxson et al., 
1979). Additional antisera were prepared to the purified albumins of 
P. chlorobryonis from McCormick County, South Carolina, P. albagula 
from Garland County, Arkansas, and P. jordani from Knob Mountain, 
Tazewell County, Virginia (reported in Maha et al., 1983). The albumins 
used to raise antisera were obtained from salamanders known to be 
homozygous for alleles 6 (P. chlorobryonis), d (P. glutinosus), and e (P. alba- 
gula) from populations 21, 70, and 126, respectively (Appendix 1). All 
antisera were produced in New Zealand White rabbits over a 13 week 
period according to procedures outlined by Maxson et al., (1979). All 
antisera were judged to be directed solely to serum albumin as evi- 
denced by a single precipitin arc in immunoelectrophoresis when tested 
with salamander plasma. Additionally, results of MC’F tests with 
purified albumin and whole plasma were indistinguishable. 
Plasma from individual salamanders from 135 populations of the 
P. glutinosus complex was compared in MC’F tests with the three 
antisera according to procedures described by Champion et al., (1974). 
For populations where gel electrophoresis had detected albumin poly- 
morphism, only individuals homozygous for each allele were used in 
MCF tests. Heterozygotes were not used in this study. When the same 
individual salamander was compared to an antiserum several times, the 
immunological distance (ID) measured was within +2 units of an 
average ID. This experimental variation in ID is the same order of 
magnitude found in previous studies and reflects the lower limit of 
resolution of the MC’F methodology (Maxson and Maxson, 1979, 1986). 
The 135 populations studied are the same as described in Part 1 and 
voucher specimens from most localities will be deposited in the collec- 
tion of the National Museum of Natural History. Nei distance (Nei, 
1972) estimates of genetic differentiation are from Part 1. The phy- 
logenetic tree in Figure 18 was constructed using a modification 
(Maxson et al., 1979) of the distance Wagner procedure (Farris, 1972). 


RESULTS AND DISCUSSION 


The 21 hour titer for the four antisera varied from 2,000 to 4,100, with 
an average of 3,000. The slope of all four antisera was 400. Both 
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parameters are typical of titers and slopes reported in other amphibian 
studies (Maxson et al., 1979; Maxson, 1984). 


P. glutinosus 
P. albagula 
P. chlorobryonis 


P. jordani 


3 2 ] 0 
X ID 


Fig. 18. Phylogeny depicting branching relationships among species of the 
Plethodon glutinosus complex for which albumin antisera were available, and P. 
jordan. Details of tree construction are in the text. Average immunological 
distances between branches of tree are indicated. 


Table. 7. Immunological comparisons (given in immunological distance units) 
of three samples of the Plethodon glutinosus complex and one sample of P. jordani 
(from Maha, et al., 1983). 


Antisera 
P. glutinosus complex P. jordani 
Antigen NY AR SC VA 
P. glutinosus complex 
NY (P. glutinosus) 0 1 4 5 
AR (P. albagula) 1 0 6 6 
SC (P. chlorobryonis) 2) 2 0 4 
P. jordan 
VA 5 9 8 0 


The reciprocal data from comparisons of the four antisera are pre- 
sented in Table 7. The average deviation from perfect reciprocity is 2 
units. The percent standard deviation from reciprocity (Maxson and 
Wilson, 1975) normally calculated for MC’F studies is 26.5 percent. 
This value is atypically high compared to earlier studies where these 
values typically average 10-15 percent. The high 26.5 percent is an 
artifact of the low ID values measured between these Plethodon popula- 
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tions. A similar situation was reported in a study of albumin evolution 
in populations of marsupial tree frogs (Scanlan et al., 1980). Because 
the average sensitivity of immunological distance measurements is +2 
ID (Maxson and Maxson, 1979), the nonreciprocity of this matrix 
becomes trivial. 

Average IDs were computed from the reciprocal values presented in 
Table 7. These values were used to draw the tree shown in Figure 18. 
The goodness-of-fit statistics for the tree are: standard deviation (Fitch 
and Margoliash, 1967) = 10.1 percent; f (Farris, 1972) = 2.0 percent; F 
(Prager and Wilson, 1976) = 7.4 percent; cophenetic correlation coeffi- 
cient (Sokal and Rohlf, 1962) = 0.959. All of these parameters indicate 
a reasonable fit of the data to the tree. For comparison, typical Fitch- 
Margoliash standard deviations in amphibian studies average 15 per- 
cent (Maxson, 1981). 

Inspection of Table 7 indicates that the P. glutinosus and P. albagula 
populations have the most similar albumins (1 ID), and that P. jordani 
has an albumin most differentiated from the three other antisera (an 
average ID of 6). P. jordani was initially included to serve as an outgroup 
for analysis of relationships within the P. glutinosus complex. Although 
it served to root the cladogram in Figure 18, Maha et al. (1983) showed 
that the albumin of P. jordani is usually more similar to the albumins of 
some members of the P. glutinosus complex than it is to the albumins of 
P. kentucki. Because P. jordani could no longer serve as an outgroup to 
the entire P. glutinosus complex, no further comparisons were made 
with this antiserum. 

Figure 18 may be somewhat misleading if we consider more extensive 
one-way albumin comparisons between the three alleles used to prepare 
antisera (summarized in Table 8). Using antisera to P. glutinosus allele d, 
the average ID to 13 samples of P. albagula allele e is 2.6 + 0.4, whereas 
the P. albagula antisera run with 36 samples of P. glutinosus allele d gives 
an average ID of 1.6 +0.2. The grand mean is 1.9 + 0.2, in contrast to 
the value of 1 obtained between the two populations used to produce the 
antisera. Similar comparisons of the P. chlorobryonis antisera to the 36 
samples of P. glutinosus allele d yield an average ID of 4.3 +0.4, 
agreeing with the single population comparison of 4. Another alterna- 
tive that cannot be ruled out is that all three species are roughly 
equidistant from one another and would arise at a trifurcation almost 
immediately after the lineage to P. jordani diverged. It is also possible 
that all four lineages are roughly equally divergent from one another. 
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Table 8. Summary of immunological distances to three antibodies of the P. glutinosus 
complex. 


New York Arkansas South Carolina 
Species Allele N antibody (d) antibody (e) antibody (b) 
Range x+S.E. Range x+S.E. Range” x=S.E. 


P. kentucki a 2 7-8 7.5 +0.5 7-9 8.0 +1.0 8 8.0 — 
P. kentucki b 2,13 13.0 — 13 13.0 — 12-13 12.5+0.5 
P. teyahalee (3 6 1-4 28 +0.5 1-7 3.8 +0.8 ~ 4-10 7.0210 
P. glutinosus c 3 3-5 4.0 +06 3-5 4.0 +06 3-9 6.0+ 1.7 
P. kisatchie Cc 2 1-4 2 5iee 5) 13 20+1.0 3-7 5.0+2.0 
P. chattahoochee d 4 3-6 45 +06 3-4 3.5 +0.3 2-5 3.8+0.6 
P. chlorobryonis d 12 0-9 2.9 +0.7 0-9 3.2 +0.7 0-7 2.0+0.6 
P. variolatus d 7 1-7 a2 =0:8. 2-5 2.77+0.4 0-6 3.0+0.8 
P. glutinosus d 36 «0-4 0.41+0.2 04 16 +02 0-10 4.3+0.4 
P. massissippr d ll 0-5 20 =O oes 2.5°20.2 Psat 4 a= 0.7 
P. cylindraceus d 3 1-3 25 EQIP e235 2.7'+03 "FR 60= 1.0 
P. kiamichi e 2 2-4 SO = LON 0-2 10+10 68 7.0+1.0 
P. grobmam e 10 2-5 3.5 +0.3 0-3 12;+0.3 412 7.3+0.8 
P. cylindraceus e 9 0-6 2.3. +0.7 0-6 3.0+06 413 8.1+0.9 
P. albagula? e 13 0-5 26+04 02 0.31+0.2 414 8.1+0.8 
P. aureolus e 2 7 2.0 — 0-2 1020.7 2 1.5+0.5 
P. ocmulgee h 4 6-7 e=20!3" "S=5 4.5 +0.5" ' GY "S32 1.1 


1 Sample size one less because homologous comparison omitted. 
2 Includes samples 92 and 116-127. 


The data presented in Table 7 show that the P. chlorobryonis antibody 
usually gives higher ID values in reciprocal tests. The data from 
comparisons of all populations to the three antisera are given in 
Appendix 3. Table 9 presents average IDs for individual alleles within 
each species. The same pattern of higher values with P. chlorobryonis is 
observed. Thus, ID values involving comparisons to this antibody may 
be slightly inflated. 

Individual salamanders from 135 populations of slimy salamanders 
were comipared to our reference panel of three antisera (Appendix 3). 
In populations where there were multiple albumin alleles, only individ- 
uals homozygous for each different allele were sampled. The IDs for all 
individuals tested are summarized in Table 9, along with the designa- 
tions of the electrophoretically identified alleles. When this study was 
initiated, all populations were considered a single species, P. glutinosus. 
Based on the analysis in Part 1, 16 species are now recognized. 
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Table 9. Average immunological distances of 135 populations of the 
P. glutinosus complex to reference antisera. 


P. glutinosus P. albagula P. chlorobryonis 


(alleles) (d) (e) (b) 
P. teyahalee 
= 6 (oc) 2.8 3.8 7.0 
P. chattahoochee 
= 4 (d) pti) 3.5 3.8 
N= 1 (e) 6.0 4.0 20) 
P. chlorobryonis 
N= 1 (b) 3.0 2.0 0.0 
N=12 (d) 29 ae 2.0 
N= 1 (g) 9.0 10.0 12.0 
P. variolatus 
N= 7 (d) 3,4 2.7 3.0 
P. ocmulgee 
N= 1 (d) 3.0 0.0 7.0 
N= 4 (h) 6.5 4.5 8.3 
P. glutinosus 
= By kb) 5.0 4.0 4.0 
N= 3 (©) 4.0 4.0 6.0 
N = 36 (d) 0.4 1.6 4.3 
N= 1 (e) 1.0 1.0 4.0 
N= 2. (h) 7.0 5.0 8.5 
P. kiamichi 
N= 2 (e) 3.0 1.0 7.0 
P. mississippi 
N=11 = (d) oF 25 8.5 
P. kisatchie 
N= 2 (©) Bell 2.0 5.0 
N= 1 (d) 5.0 3.0 7.0 
P. sequoyah 


N= 1 (d) 2.0 3.0 14.0 
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P. glutinosus P. albagula P. chlorobryonis 
(alleles) (d) (e) (b) 


P. grobmani 

N=10 (e) 3.5 1.2 7.3 

N= 1 (g) 7 6.0 8.0 
P. cylindraceus 

N= 1 (©) 4.0 5.0 10.0 

N= 3 dd) 2.3 P| 8.0 

N= 9 (e) 2.3 3.0 8.1 

N= 1. (h) 7.0 6.0 11.0 
P. albagula 

IN =. RO“) 5.0 0.0 6.0 

N=13 (e) 2.6 0.3 8.1 

N= 2S1G) 4.0 5.0 13.0 
P. savannah 

N= keh) 7.0 4.0 5.0 
P. aureolus 

N= 1 (@) 3.0 2.0 0.0 

N= 2 (e) 2.0 1.0 1.5 
P. kentuchi 

N= 2 (a) 7.5 8.0 8.0 

N= 2 (b) 13.0 13.0 12.5 


The highest immunological distances are found in the comparisons of 
P. kentucki allele b with all three antibodies. The other common P. ken- 
tucki allele a is also among the most different from the panel of 
antibodies. This result agrees with the finding that P. kentucki is 
genetically the most divergent species of the P. glutimosus complex (Part 1). 

The P. chlorobryonis antibody was prepared from animals homozygous 
for allele b. This allele is absent in most other populations of P. chlo- 
robryonis, but its mean ID to the common allele d of P. chlorobryonis (2.0) 
is lower than to most other alleles in other species. With the exception of 
P. aureolus (discussed below), only the comparisons to allele d in P. 
chattahoochee and P. variolatus are less than 4.0. The latter two species are 
genetically most closely related to P. chlorobryonis (Part 1). 

Throughout its range, P kentucki is sympatric with P. glutinosus. 
Despite the difficulty in distinguishing the two species using mor- 
phological criteria, P. kentucki has two common albumin alleles that are 
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distinct from that of P. glutinosus by both electrophoretic and immu- 
nological criteria (Maha et al., 1983, Highton and MacGregor, 1983). 
Throughout its range, P. aureolus is sympatric with P. teyahalee and the 
two species are sympatric with P. glutinosus at one locality in Polk 
County, Tennessee (Highton, 1984). All other species replace each other 
geographically (Part 1). 

Surprisingly, the immunological distances of the albumins of all three 
populations of P. aureolus to the panel of antibodies are low (0—3 ID). 
This result is unexpected in light of the considerable amount of 
electrophoretic genetic differentiation of this species compared to all 
the others (Part 1). Perhaps the P. aureolus albumin incurred signifi- 
cantly fewer substitutions than some of the other loci screened electro- 
phoretically. 

Tables 8 and 9 show considerable albumin variation within several of 
the species of the P. glutinosus complex. Within species, there are alleles 
that differ (by immunological criteria) from those of the reference 
antisera. For example, among the 40 populations of P. glutinosus, five 
albumin electromorphs are observed. Alleles designated b and c have 
essentially the same pattern of immunological cross-reactivity as do 
alleles d and e. Allele A exhibits still a third pattern. We interpret this 
amount of differentiation as indicative of a considerable amount of 
evolutionary divergence. Such albumin differentiation has also been 
found in other amphibians with little or no discernable morphological 
variation (Maxson et al., 1977; Maxson, 1981; Maxson, unpub.). 

Twenty-four of the 29 populations of P. glutinosus from the northern 
portion of its range (NY, NJ, PA, MD, WV, VA, OH, IN, IL, and KY) 
have immunological distances of 0 to the homologous New York 
antibody. These northern populations are all genetically closely related. 
This species very likely extended its range northward from south- 
eastern Kentucky since the last glaciation (Part 1). These results might 
be expected if the source populations were monomorphic for a single 
albumin. The remaining populations of P. glutinosus with the 6, c, e, and 
h alleles occur in the southern part of the range of the species (Figs. 1 & 
11). The Nei genetic distances (Appendix 3) between the New York 
population and these southern populations are higher (D = 0.09) than 
the Nei genetic distances to the northern populations (D = 0.03). 

Another striking feature of these data is that alleles with the same 
electrophoretic mobility (same letter designation) do not always behave 
identically to our panel of antisera. For example, allele d from New York 
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P. glutinosus is immunologically different from allele d found in P. chatta- 
hoochee, P. chlorobryonis, P. variolatus, P. ocmulgee, P. mississippi, P. kisatchie, 
P. sequoyah, P. cylindraceus, and P. aureolus. Without additional antisera to 
these other species, it is not possible to determine if all d alleles are 
identical to one another. However, by examining their pattern of cross 
reactivity to the other two antisera (compare the d allele of P. chatta- 
hoochee, P. chlorobryonis, and P. variolatus with the d allele of P. mississippi 
and P. cylindraceus, Tables 8 & 9) it appears that there are at least two 
additional “d” alleles. These results are not surprising in light of the 
genetic distances among these species. 

Another example involves allele e for which the antiserum to P. alba- 
gula was prepared. The e alleles of P. grobmani and P. kiamichi show a 
common pattern of cross reactivity to P. albagula. However, allele e in 
P. aureolus shows yet another pattern of cross reactivity. 

It would be expected that some mutations might not change electro- 
phoretic mobility (Reichlin et al., 1966; Cocks and Wilson, 1969; Prager 
and Wilson, 1976; King and Wilson, 1975). Sequential gel electro- 
phoresis (Coyne, 1982) might also have detected some of these addition- 
al alleles. Maiorana (unpublished data) found that the mobility of allele 
din P. glutinosus was different from the d allele in P. chlorobryonis when 
compared using the method of Smithies (1959). King and Wilson (1975) 
estimated that only 28 percent of all amino acid substitutions that occur 
are resolved by standard gel electrophoresis but Coyne (1982) suggests a 
higher estimate (about 50 percent). The fact that banding patterns seen 
on electrophoretic gels are phenotypes, and thus do not necessarily 
reflect the complete underlying genotype, has been discussed by others 
(Allendorf, 1977; King and Ohta, 1975). It has been demonstrated, 
however, that MC’F can detect those substitutions that do not change the 
overall net charge of the protein (Cocks and Wilson, 1969; Ibrahimi et 
al., 1980). Our results also show the ability of MC’F to detect additional 
alleles. 

Not all alleles show this confusing pattern. For example, allele h in 
P. ocmulgee (populations 31-34), P. glutinosus (populations 72, 73), P. cyl- 
indraceus (population 111), and P. savannah (population 128) all show 
similar patterns of cross reactivity to the panel of antisera. It would be 
interesting to find out whether allele h is the same in these genetically 
rather different but geographically parapatric species. 

Some alleles may not be distinguished by immunological analysis. For 
example, MC’F did not distinguish between alleles d and e in P. glutinosus 
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or dand ¢ of P. cylindraceus (Tables 8 & 9). This inability to discriminate 
between alleles is not surprising. An antiserum to each allele in the 
study as well as many more populations would probably be needed to 
discriminate between all of the different albumins. 


Jordanl 
(Knob Mt) 


glutinosus \c) 


cylindraceus (d) 


chlorobryonis (d) cylindraceus (e) 
variolatus (d)" 
mississippi (d) 
glutinosys (d) 
aureolus (e) 


fae) 

) 
kiamichi (e) 
grobmani (e) 


kentucki (a) 


5 


ocmulgee (h) 


Fig. 19. Network of inferred relationships among species in the P. glutinosus 
complex, indicating all of the common alleles detected by gel electrophoresis. 
The distances along the branches are immunological distance units. 


An analysis of the correlation of immunological distance and Nei 
genetic distance was performed because there are 431 comparisons. The 
correlation coefficient was 0.585, considerably lower than that observed 
in earlier studies where correlations of between 0.8 and 0.9 were 
reported (Maxson and Maxson, 1979; Highton and Larson, 1979; Wyles 
and Gorman, 1980). However, all earlier studies involved IDs ranging 
from 0-80 ID, with very few IDs less than 10 units. This study, on the 
other hand, involved IDs all less than 14. Despite the small range of 
immunological distance in these comparisons, the correlation was still 
significant at the 0.01 level. The slope of the linear regression line for 
these data is 15.5, also somewhat lower than that in previous studies of 
salamanders: Hydromantes (28), Pseudoeurycea (23), and Plethodon (24) 
(Maxson and Maxson, 1979). This variation in slope should serve as a 
cautionary warning to investigators who use Nei genetic distances to 
estimate time of divergence from albumin immunological studies. The 
apparent variation in the slope appears to be attributable to the species 
studied as well as to the suite of loci compared. 

Figure 19 is a hypothetical network showing possible relationships 
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among the more common albumin alleles. This network is reasonably 
consistent with the MC’F data. Because we had antibodies to only three 
alleles, many of the hypothetical relationships are tentative. In general, 
the relationships agree with those described in Part 1. Examples are the 
similarity of P. chattahoochee, P. chlorobryonis, and P. variolatus, the isolated 
position of P. kentucki, and the approximately equidistant relationship of 
the remaining species. Amino acid sequencing of the various albumins 
of the P. glutinosus complex would be needed to better understand the 
evolutionary relationships of these proteins. 


SUMMARY 


This immunological analysis of albumin differentiation in the P. glu- 
tinosus complex has shown that quantitative MC’F can detect a consider- 
able amount of the variation present in closely related species. As 
suggested earlier (Maxson and Maxson, 1979) electrophoretic analysis is 
the methodology of choice for studying closely related populations. 
Although MC’F can detect single amino acid substitutions, even in the 
absence of a charge change, an antiserum to each genetically differenti- 
ated population is needed for more definitive work. With the three 
antisera used in this study, we could identify populations as different 
but could not characterize taxonomic groupings as is possible with 
electrophoretic analysis (Part 1). By combining MC’F and electrophore- 
tic analyses of salamander albumins, we have demonstrated more 
underlying genetic variation than is resolvable by either technique 
alone. 
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Appendix 3. Immunological distances of the Plethodon 


glutinosus complex 


Population 
County 
or 
No. State parish 
P. teyahalee 
Ne CN Monroe 
2 TN Monroe 
3 NC Macon 
4 NC Haywood 
5 NC Graham 
6 NC Madison 


P. chattahoochee 
7 GA Union 
8 GA Fannin 
9 GA Towns 
10 GA Rabun 


P. chlorobryonis 
11 SC Chesterfield 


12 SC Florence 

13 NC Columbus 

14 NC Craven 

15 NC Halifax 

16 VA Southampton 

17 VA Dinwiddie- 
Sussex 

18 NC Tyrrell 

19 NC Pender 

20 GA Habersham 

21 SC McCormick 


22 SC Aiken 


Allele 


aaaadao a (ot oP oF ob oF @) i (a) (a @ @ 


aaraa ea 


New York 


(pop. 70) 


D 


ID 


Coe BPO Oo PR OO PP a HD Ov P 09 Oo ee BB OO P NO 


— © OOO & 09 


Antibody 

Arkansas 

(pop. 126) 

D ID 
0.37 5 
0.30 + 
0.35 1 
0.34 di 
0.32 3 
0.31 3 
0.27 3 
0.30 3 
0.30 4 
4 
0.24 4 
5 

0.25 10 
0.27 ] 
0.27 3 
0.29 1 
0.28 4 
0.28 4 
0.28 4 
0.28 1 
0.27 3 
0.29 9 
2 
0.31 0 
0.30 3 


South Carolina 


(pop. 21) 
D ID 
0.49 10 
0.42 4 
0.47 6 
0.47 5 
0.44 7 
0.43 10 
0.29 4 
0.31 4 
0.20 5 
2 

0.19 9 
1 

0.17 12 
0.14 2 
0.15 0 
0.15 7 
0.16 ] 
0.16 1 
0.15 a 
0.15 1 
0.14 3 
0.09 5 
0 

0.06 0 
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Population Allele Antibody 
County New York Arkansas South Carolina 
or (pop. 70) (pop. 126) (pop. 21) 
No. State parish D ID D ID D ID 
P. variolatus 
23 SC Charleston- d 0.35 2 0.30 3 0.19 4 
Dorchester 
24 SC Jasper d 0:35," 2 0.45: 5 0.19 ] 
25 SC Allendale d 0.36 2 0.33 2 0.20 0 
26 SC Jasper d O34. 47°" O30) a2 0.24 3 
27)SG Berkeley d 0.26 4 0.24 2 0.23 4 
28 GA Chatham d 0.18 a 0.14 3 0.20 6 
29 GA _ Effingham d 0.16 1 0.24 2 0.28 3 
P. ocmulgee 
30 GA _ Bacon d 0.20 3 0.27 0 0.36 7 
31 GA Bulloch h 0.17 a 0.20 5 0.42 11 
32 GA _—~Wheeler h 0.19 6 0.23 5 0.38 9 
33 GA Upson h 0.205 7 0:25 “5 kaa | 
34 GA __— Long h 0.26 6 0.30 3 0.39 6 
P. glutinosus 
35 IN Moore c 0.14 5 0.17 4 0.39 6 
36 AL Etowah b 0.11 5 0.15 4 0.99 4 
c 4 3 3 
37 TN _ Bledsoe c 0.11 3 0.12 5 0.36 9 
38 TN Polk d 0.09 0 0.16 0 0.36 3 
39 KY Muhlenberg d 0.04 1 0.20 2 0.42 1 
40 IN Crawford d 0.05 0 0.19 2 0.42 6 
41 IL Union d 0.03 0 0.21 3 0.37 a 
42 IL Pope d 0.04 0 0.21 3 0.38 5 
43 KY Warren d 0.05 0 ON? 2 0.37 0 
44 IN Parke d 0.03 0 0.16 1 0.41 6 
45 TN _ Blount d 0.03 0 0.16 ] 0.37 5 
46 AL Cleburne d 0.06 1 0.13 0 0.32 5 
47 KY McCreary d 0.03 0 0.15 0 0.34 a 
48 NJ Union d 0.03 0 0.16 3 0.37 6 
49 VA Giles d 0.04 0 0.15 2 0.32 2 
50 OH Meigs d 0.04 0 0.15 3 0.31 5 
51 VA Russell d 0.05 0 0.15 1 0.32 7 
52 IN Jackson d OOS, “0 4 O15) 3 0.36 8 
53 KY Harlan d 0.03 0 0.16 0 0.31 2 
54 OH ~ Brown d 0.02 0 0.15 2 0.32 a 
55 MD _ Frederick d 0.02 0 0.16 1 0.35 7 
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Population Allele Antibody 
County New York Arkansas South Carolina 
or (pop. 70) (pop. 126) (pop. 21) 
No. State parish D ID D ID D ID 
56 PA Cambria d 0.02 0 0.16 1 0.33 3 
57 WV _ Preston d 0.02 1 0.16 2 0.33 5 
58 PA Susquehanna d 0.01 Or VOTS. 2 0.33 3 
59 NY ‘Tompkins d 0.04 2 0.15 l 0.30 4 
60 PA Beaver d 0.03 2 0.16 4 0.32 2 
61 PA Snyder d 0.03 0 0.16 3 0.31 4 
62 PA Cameron d 0.03 0 0.16 | 0.31 4 
63 KY Pike d 0.03 1 0.15 1 0.31 3 
64 WV Clay d 0.03 0 0.15 1 0.31 6 
65 VA Wise d 0.03 0 0.15 1 0.31 3 
66 KY Harlan d 0.03 0 0.16 2 0.30 3 
67 OH Wayne d 0.05 0 0.16 3 0.30 4 
68 KY Scott d 0.03 0 0.18 0 0.33 6 
69 PA Bedford d 0.01 0 0.17 2 0.34 7 
70 NY Ulster d - — 0.18 1 0.38 4 
71 TN Montgomery d O07 4 0 ~ O18 0.35 a 
723GA / ‘Henry d 3 3 1 
h 0.09 7 0.16 5 0.25 7 
13° AL Macon d 4 1 10 
h 0.11 - 0.14 P 0.33 10 
74 TN — Cocke e 0.10 1 0.15 ] 0.44 4 
P. kiamichi 
75 OK Le Flore e 0.30 4 0.16 0 0.36 6 
76 OK Le Flore e 0.30 2 0.17 2 0.33 8 
P. mississippi 
TAL Tuscaloosa d 0.16 0 0.14 3 0.29 11 
78 AL Lawrence d 0.19 5 Only, 1 0.35 10 
79 MS Tishomingo d 0.19 1 0.16 3 0.33 8 
80 TN Henderson d 0.16 1 0.26 1 0.45 5 
81 TN Shelby d 0.16 3 O27. 2 0.45 9 
82 MS Winston d 0.21 3 0.21 2 0.36 8 
83 MS Scott d 0.22 3 0.21 3 0.38 10 
84 MS Lowndes d 0.19 4 0.16 3 0.35 10 
85 LA Washington d 0.22 4 0.20 3 0.38 11 
86 MS Forrest d 0.23 4 0.22 2 0.39 8 
87 AL Clarke d 0.27 2 0.25 2 0:37 4 


P. hisatchie 
88 AR Union d 0.33 5 0.25 3 0.40 7 
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Population Allele Antibody 
County New York Arkansas South Carolina 
or (pop. 70) (pop. 126) (pop. 21) 
No. State parish D ID D ID D ID 
SOnEAy = (Winn c 0405" 4 80:25 1 0.43 3 
90 LA Grant c 0.40 OD. 23 0.43 7 
P. sequoyah 
91 OK = McCurtain d 0.29.5 2 .0.20" 43 0.40 14 
P. albagula 
92 AR Montgomery c 0.21 5 0.09 0 0.29 6 
e 3 0 6 
P. grobmani 
93 FL Hillsborough =e 0.95 3 0.15 3 0.36 4 
g 7 6 8 
94 FL Columbia e 022993 010 2 0.34 8 
95 GA Lanier e O2atar 5 OL: 20 0.34 8 
96 FL Lake e 0265" 4 0:09 | 2 0.37 7 
97 AL Butler e O14, 4 0:12 1 0.38 12 
98 FL Leon e 0.21 4. Ori 1 0.33 5 
99 GA __ Schley e 0.21 3 = 009) 0 0335 6 
100 AL Barbour e 0.22 53 © 0038 0 0.29 9 
101GA Terrell e 0.21 2 = 0:09 ] 0.32 5 
102FL ‘Jackson e 0:25 (4° Ol 2 0.32 9 
P. cylindraceus 
103NC Montgomery e O28 0 O15 3 0.29 11 
104 VA Grayson- e O29" A ORG Fees 0.33 q 
Wythe 
105 VA Fairfax e O38 ~2. “0:22 70 0.42 = 
106 VA Bedford e 036 "0 | O21 2 0.40 8 
107 VA Botetourt & Oba =3 ~10:20° PSs 0.39 8 
108 VA Augusta e 36°" 0 ~ 020 3 0.40 6 
109WV Hampshire e 0.33 oS - OSS. 0.38 = 
110 VA Pittsylvania e O33a7 . 3 ONS 2 0.37 5 
111 NC Burke e 6 6 13 
h 0.32 7 0.19 6 0.34 1 
112SC Chester d 0.21 I, 08S « 23 0.30 10 
113.NC- = Madison- d Oa oo | O22) eo 0.31 7 
™N Unicoi 
114NC- — Mitchell- d O23" 3: "O22 2 0.31 7 
T™N Unicoi 


ae 


115 NC Henderson Cc 0.32 0.24 5 0.34 10 
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Population 
County 
or 
No. State parish 
P. albagula 
116TX Hays 
117AR Pope 
118MO _ Iron 
119MO  +Wayne 
120MO ‘Miller 
121 AR Stone 
m2 Tx ‘Travis 
1231 X Kerr 
124OK  Sequoyah 
125 OK Adair 
126 AR Garland 
127 AR Dallas 
P. savannah 
128GA Richmond 
P. aureolus 
129 TN Monroe 
130 TN Monroe 
131TN Polk 
P. kentucki 
132 VA Wise 
133 KY Harlan 
134 KY Pike 
135 KY Harlan 


Allele 


oon n rnrnrnarrm md O® =: 


vnrao Te 


New York 


(pop. 70) 


D 


0.35 


0.26 
0.32 
0.32 
0.30 
0.29 


ID 


— e OF eS OF O NP PLP BP OO #P OO 


Antibody 
Arkansas 
(pop. 126) 
D ID 
] 
0.20 5 
0:07 ~2 
0.12 od 
012-— 40 
Ont 0 
0.11 0 
0. 144/980 
Cis “0 
0.05 O 
O07, 0 
0.01 0 
0.21 4 
0.24 O 
O27 a2 
0.34 2 
0.43 
0.42 3 
0.41 7 
0:39: -"g 
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South Carolina 


(pop. 21) 
D ID 
5 

0.31 13 
0.37 9 
0.38 1] 
0.38 14 
0.37 1 
0.38 4 
0.33 7 
0.38 8 
0.37 6 
0.39 1] 
0.31 6 
0.32 i 
0.36 GB 
0.50 2 
0.57 ] 
0.56 0 
0.51 12 
0.54 ills} 
0.53 8 
0.53 8 
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81 
Polymorphism, 82 
Allospecies, 52 
Alpha-Glycerophosphate dehydrogenase 
(a-Gpd), 6, 39 
Amino acid sequence difference, 81, 82 
Antibody, albumin, 82-92 


Batrachoseps, 50 
Bolitoglossa occidentalis, 48 
Bolitoglossa rufescens, 48 
Buffers, 6 
Lithium hydroxide, 6 
Poulik, 6 
Tris-citrate, pH 6.7, 6 
Tris-citrate, pH 8.0, 6 
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Tris versene borate, 6 


Computer program 
BIOSYS-1, 7, 23, 31 
Taxan, 7 
Contact zones, 52, 53 
Cophenetic correlation coefficient, 22, 23 
Cylindrosoma glutinosum, 62 


Distance 
Build-up of genetic, 12, 50 
Cavalli-Sforza and Edwards, 23 
Edwards, 23 
Immunological, 83-86, 91, 145-149 
Nei standard genetic (D), 6, 12, 13, 23, 

26—30, 51, 83, 91, 145-149 

Prevosti, 23 


Rogers, 23 
Drosophila, 47 
Drosophila willistont group, 47 


Electrophoresis, 82 
Acrylamide disc, 6 
Sequential, 90 
Starch gel, 6 

Esterase (Est), 6, 9, 39 

Evolution 
Differential rates, 39, 42 
Morphological, 12 
Protein, 12 


Fumarase (Fum), 6 


Gel, 6 
Acrylamide disc, 6 
Starch, 6 
Genetic 
Data, 4, 105-144 
Differentiation 
Species level, 13 
Structural gene, 49 
Distance. See Distance: Nei standard 
genetic 
Identity (Nei’s J), 6-8 
Loci, substitutions within, 39-43 
Variability, 43-45 
Glutamate dehydrogenase (Gdh), 6, 9, 13, 
42 
Glutamic oxaloacetic transaminase (Got), 
6, 39 
Grouping of populations, 9-23 
Groups. See under species names 
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Hennigian methods, 49 
Heterozygosity (H), 7, 43-45 
Hydromantes, 91 
Slope of linear regression between JD 
and D, 91 


Interpopulation relationships, 31-37 
Pleistocene effects on, 37-39 
Isocitrate dehydrogenase (Idh), 6, 9 


Lactate dehydrogenase (Ldh), 6, 39 
Leucine aminopeptidase (Lap), 6 
Localities, list of, 101-104 


Malate dehydrogenase (Mdh), 6, 9, 39 
Micro-complement fixation (MC’F), 81, 
82, 90, 92 

Molecular clock, 12, 26 
Molecular evolution, 82 

Equal rates of, 14, 22 
Monomorphic locus, 4, 6 
Morphological evolution 

Convergent, 82 

Relation to molecular evolution, 82 


National Museum of Natural History 
(USNM), 9 
Null alleles, 9 


Peptidase (Pep), 6, 9 


Phenogram, UPGMA, 7-9, 14, 22, 23, 25, 


a3, 25, 30, 43 
Phosphoglucomutase (Pgm), 6 
Phosphoglucose isomerase (Pgi), 6 
Phylogeny, errors in reconstruction, 20 
Pleistocene effects on interpopulation 

relationships, 26, 34, 37-39 
Plethodon 

Previous genetic variation studies, 4 

Life history study, 47 

Slope of linear regression between JD 
and D, 91 

Species groups, 12 

Plethodon albagula (group 13), 18, 19, 35, 

36, 38, 46, 71-73, 78 

Albumin antibody, 83, 84, 87, 88 

Genic variation, 131-134, 137-140 


Heterozygosity estimates, 43, 134, 140 


Immunological distances, 84—86, 88, 
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91, 145-149 
Plethodon aureolus (group 15), 3, 9, 14, 20, 
25, 26, 37, 45, 46, 50, 52, 53, 75 
Genic variation, 141-144 
Heterozygosity estimates, 43, 144 
Immunological distances, 86, 88-91, 
149 
Plethodon caddoensis, 39, 43, 51 
Plethodon chattahoochee (group 2), 15, 33, 
34, 46, 55-58 
Genic variation, 105—109 
Heterozygosity estimates, 43, 109 
Immunological distances to, 86, 87, 90— 
42,145 
Plethodon chlorobryonis (group 3), 15, 16, 
34, 36, 37, 46, 58, 59 
Albumin antibody, 83-85, 87, 88 
Genic variation, 109-112 
Heterozygosity estimates, 43, 112 
Immunological distances, 84—87, 90— 
92, 145-149 
Plethodon cinereus, 31, 34, 37 
Plethodon cylindraceus (group 12), 16, 18, 
19, 33, 35, 36, 38, 46, 70-72 
Genic variation, 134-137 
Heterozygosity estimates, 43, 137 
Immunological distances, 86, 88, 90, 
91, 148 
Plethodon fourchensis, 39, 51 
Plethodon glutinosum, 61 
Plethodon glutinosus (group 6), 3, 9, 14, 17, 
34-36, 38, 46, 52, 54, 61-63, 75, 76 
Albumin antibody, 82-84, 87, 88 
Genic variation, 115—127 
Heterozygosity estimates, 43, 118, 121, 
1245127 
Immunological distances, 84-87, 89— 
91, 145-149 
Species complex, 3, 7, 10, 12, 15, 16, 31, 
42, 43, 45, 50, 51, 53, 54, 76, 82, 84, 
92 
Species group, 12, 39, 42, 43, 45 
Plethodon glutinosus albagula, 71. See also 
Plethodon albagula 
Plethodon glutinosus chlorobryonis, 58. See 
also Plethodon chlorobryonis 
Plethodon glutinosus glutinosus, 54, 55, 62. 
See also Plethodon glutinosus 
Plethodon glutinosus grobmani, 69. See also 
Plethodon grobmani 
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Plethodon grobmani (group 11), 18, 35, 36, 
39, 46, 69, 70 
Genic variation, 131—134 
Heterozygosity estimates, 43, 131-134 
Immunological distances, 86, 88, 90, 
91, 148 
Plethodon jordani, 14, 15, 26, 31, 33, 39, 51, 
54, 55 
Albumin antibody, 83-86 
Complex, 7 
Immunological distances, 84, 85, 91 
Plethodon jordani kentucki, 75 
Plethodon jordani rabunensts, 15, 55 
Plethodon jordan: teyahalee, 54 
Plethodon kentucki (group 16), 3, 9, 20, 23, 
25, 26, 37, 45, 46, 50, 52, 53, 62, 75, 
76 
Genic variation, 141-144 
Heterozygosity estimates, 43, 144 
Immunological distances, 85, 86, 88, 
91, 92, 149 
Plethodon kiamichi (group 7), 17, 35, 36, 
46, 62, 64, 65 
Genic variation, 127-131 
Heterozygosity estimates, 43, 131 
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91, 147 
Plethodon kisatchie (group 9), 18, 35, 36, 
38, 46, 67, 68 
Genic variation, 131-134 
Heterozygosity estimates, 43, 134 
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91, 147, 148 
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36, 38, 39, 46, 65-67 
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91, 147 
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35, 39, 42, 46, 60, 61 
Genic variation, 112—115 
Heterozygosity estimates, 43, 115 
Immunological distances, 86, 87, 90, 
91, 146 
Plethodon ouachitae, 17, 35, 39, 43, 51 
Plethodon petraeus, 42 
Plethodon punctatus, 39, 43 
Plethodon savannah (group 14), 20, 36, 37, 
46, 73, 74 
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Genic variation, 131-134 
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Immunological distances, 87, 90, 91, 
148 
Plethodon serratus, 31 
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Plethodon shermani rabunensis, 55 
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25, 31, 33,35, 36, 38, 40, Sabana 
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Hybrid origin, 33 
Immunological distances, 86, 87, 91, 
145 
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46, 59, 60 
Genic variation, 112—115 
Heterozygosity estimates, 43, 115 
Immunological distances, 86, 87, 90, 
91, 146 
Plethodon variolosum, 59 
Plethodon wehrlei, 39, 43 
Group, 39, 43 
Plethodon yonahlossee, 39, 51 
Populations studied, 4 
Selection of, 4 
Sites, 5, 101-104 
Protein-buffer combinations, 4, 6 
Protein-1, 6 
Protein-2, 6, 9, 14, 42 
Protein-3, 4, 6, 39 
Pseudoeurycea, 91 
Slope of linear regression between JD 
and D, 91 


Rana pipiens complex, 53 
Relationships among groups, 23-37 
Reproductive isolation, 50 


Salamandra cylindracea, 16, 18, 70 
Salamandra glutinosa, 17, 61 
Salamandra melanoleuca, 62 
Salamandra variolata, 16, 18, 59 
Sample size, 4 
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Six-phosphoglutonate dehydrogenase 
(6-Pgd), 6, 39 
Size, variation among groups, 45-47 
Species, 50-53 
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85 
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and D, 91 
Subspecies, 47, 50-53 
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Substitutions within loci, 39-43 
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Taxonomy, 47—76 
Transferrin (Trf), 6, 7, 9, 39 
Tree, phylogenetic, 22, 81, 84 
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